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y-RAYS PRODUCED BY SLOW NEUTRON CAPTURE IN 
BERYLLIUM, CARBON, AND NITROGEN! 


By B. B. Kinsey, G. A. BARTHOLOMEW, AND W. H. WALKER* 


Abstract 


The y-rays produced by slow neutron capture in beryllium, carbon, and 
nitrogen have been investigated with a pair spectrometer. Only one y-ray was 
observed to follow capture in beryllium and in carbon; nine were observed from 
nitrogen. The neutron binding energies computed from the energies of these 
radiations are: Be!®, 6.797 + 0.008 Mev.; C'’, 4.948 + 0.008 Mev.; N", 10.8234 
0.012 Mev. The nitrogen y-rays can be fitted into a term scheme containing 
the previously known excited states of N“, obtained from the study of (d, p) 
reactions, together with one excited state at 9.16 Mev., hitherto unreported. 
A comparison of the intensities of the nitrogen and beryllium radiations shows 
that the radiative capture cross section of nitrogen is at most 160 millibarns. 


Introduction 


A preliminary report on the y-rays produced by thermal neutron capture in 
beryllium, carbon, and nitrogen has appeared elsewhere (15). A more detailed 
description of these measurements is presented in this paper. 


The method of measurement will be fully described in an article to be pub- 
lished shortly. The samples of material to be studied were enclosed in dural 
containers which were placed in a region of high neutron flux near the reacting 
vessel in the pile. A bismuth block was placed between the sample and the 
reactor to prevent radiation from the pile from entering the beam to be studied. 
Both the dural container and the bismuth block produce unwanted radiations 
which are superposed on the energy spectrum of the y-rays produced by the 
sample. These unwanted radiations were identified by a separate experiment 
in which an empty container was used.f 

In a survey measurement for each sample, the y-ray spectrum was explored 
from 3 Mev. to beyond that expected for the binding energy of the product 
nucleus. This measurement was followed by detailed studies of the coincidence 
peaks produced by individual radiations. The intensities of the radiations were 
such that counting times up to one hour were required at each point. The 
resolution of the spectrometer was chosen to be about 4% for the survey of the 
y-rays produced by carbon, and about 2% for all the other measurements. 


1 Manuscript received August 25, 1950. 
Contribution from the Atomic Energy ip National Research Council of Canada, 
Chalk River, Ontario, Canada. Issued as N.R.C. No. 2266. 
* Present address: Department of Chemistry, MiMaster University, Hamilton, Ontario. 
t The spectrum of the empty container (see Fig. 1) is mainly composed of aluminum radiations; 
two small bumps on either side of the ground state radiation of aluminum are due to copper and 
manganese. Magnesium does not contribute to the container spectrum in measurable amount. 
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Beryllium was irradiated in the form of beryllium oxide. Radiations from 
oxygen have not been observed, partly because the cross section for capture is 
very low—about 1.6 millibarns—and because the binding energy of the neu- 
tron in QO” is low (about 4 Mev.). 


Carbon was irradiated in the form of a solid graphite cylinder 4 in. in dia- 
meter and 6 in. long. The cylinder was cut from a graphite sample of very high 
purity. A block of chemically pure lead, 4 in. in diameter and 4 in. long, was 
placed in the experimental hole between the bismuth plug and the graphite 
sample in order to absorb the y-rays produced by neutron capture in the 
aluminum container enclosing the bismuth plug. An extra lead collimator was 
placed between the graphite and the spectrometer to ensure that no radiations 
emitted from the sides of the hole could be seen by the spectrometer. With 
these precautions aluminum radiations were prevented from entering the 
spectrometer and ensured that the only radiations below 6 Mev. which might 
be detected must be due to carbon. Other experiments have shown that lead 
emits no appreciable amount of radiation below 6.68 Mev. (16). 


For the nitrogen experiments samples of beryllium nitride, BesN2, and urea, 
CO(NH2)2, were used. Both samples were enclosed in dural containers. The 
nitrogen radiations are easily recognized from a comparison of the spectra 
produced by these compounds and by beryllium oxide. 


Beryllium 


The y-ray spectrum produced by beryllium oxide in its container and that 
of the container itself is shown in Fig. 1*, the ordinates of the latter curve being 
reduced to make the heights of the 7.7 Mev. aluminum peak the same. Within 
the statistical accuracy the two spectra coincide everywhere except for the 
broad peak at 4.17 Mev. which appears enhanced in the spectrum of the con- 
tainer, and the peak at 6.79 Mev. and its immediate vicinity. The 4.17 Mev. 
radiation is due to neutron capture in bismuth and is described elsewhere. (16). 
This radiation is not appreciably absorbed in the container and is therefore 
observed in the container spectrum with its full intersity. The transmission 
of this radiation through beryllium oxide reduces its intensity by about 40%, 
in rough agreement with the amount expected from calculation. 


Fig. 2 shows the results of a more detailed examination of the 6.79 Mev. 
radiation, the only y-ray shown in Fig. 1, which can be ascribed to beryllium. 
The energy of the y-ray is 6.794 + 0.008 Mev. When corrected for the recoil 
of the beryllium nucleus, the transition energy, which is the binding energy 
of the neutron in Be”, is found to be 6.797 + 0.008 Mev. 


This energy is in good agreement with that calculated from the energy 
balance of the Be*(d, »)Be"® reaction. For this quantity, Buechner and Strait 


* In all the coincidence spectra shown in this paper, the energy of a y-ray is given by the abscissa 
of the linearly extrapolated upper limit of its coincidence peak. The ordinates are not corrected for 
the variation of the efficiency of the spectrometer with energy. 
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Fic. 1. Survey spectrum of beryllium capture gamma-radiation, 
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(4) have found the value 4.576 + 0.012 Mev. Using the value of the disinte- 
gration energy of the deuteron obtained by Bell and Elliott (1), viz., 2.231 + 
0.007 Mev., and adding the equations: 

Be? + H? = Be” + H!+ 4.576 + 0.012 Mev. 

H' + nm = H? + 2.231 + 0.007 Mev. 


the neutron binding energy is given by:— 
Be® + mn = Be” + 6.807 + 0.014 Mev. 


Fig. 1 shows little evidence for the existence of any other radiations emitted 
by beryllium. The small bump on the low energy side of the 6.79 Mev. peak 
is substantiated by only two points, which are insufficient evidence to prove 
the existence of a y-ray. Only one bound excited state is known in Be". Its 
energy is 3.375 Mev. (4). A radiative transition from the capturing state to 
this excited state would have an energy of 3.42 Mev. With a resolution of 2%, 
the two radiations with energies of 3.375 and 3.42 Mev. would probably 
appear as a single broad peak in the coincidence spectrum. No such peak was 
observed (Fig. 1). The ratio of the efficiency for detection by the spectrometer 
at 3.4 Mev. to that at 6.8 Mev. is 0.07. If it is assumed that the 3.4 Mev. peak 
would be detectable if its intensity were 0.6 coincidence per minute above the 
background (or a total of 1.6 coincidences per minute in Fig. 1) then the ratio 
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of the intensity of the cascaded transitions to the direct transition would be 
about 0.4. Transitions to the 3.375 Mev. excited state, if they occur, must 
therefore have less than 0.4 of the intensity of the direct transition to the 
ground state. 


oa 
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Fic. 2. Detail of beryllium radiation. 


Carbon 


The result of a survey measurement of the y-rays emitted by carbon between 
3 and 5.5 Mev. is shown in Fig. 3. The 4.9 Mev. radiation is emitted in the 
direct ground state transition in C. There is some indication of a radiation 
at 5.3 Mev. It is unlikely that a coincidence peak at this point could be pro- 
duced by any of the known impurities in the graphite sample used in these 
experiments. The only impurity with a strong radiation in this region of energy 
is silicon. However, the amount of this impurity and its cross section are too 
low to account for the peak observed. No evidence for such a radiation was 
found in measurements taken with a higher resolution; however, the counting 
rates observed were so low that a coincidence peak of the height expected was 
difficult to distinguish against statistical fluctuations. 


There is some evidence in Fig. 3 for radiations with energies between 3.4 and 
3.7 Mev. These radiations must be due to impurities, for no y-rays produced 
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by the excitation of the two low-lying states of C™ (5, 17) at 3.1 and 4.0 Mev. 
can lie in this range. An upper limit may be put on the intensities of these latter 
radiations if the minimum detectable peak coincidence counting rate is as- 
sumed to be twice the probable error of the ordinates of Fig. 3 at the corres- 
ponding energy. The upper limit for the intensity of the 3.1 Mev. radiation is 
found to be about 75%, and that of the 4.0 Mev. radiation about 30% of the 
intensity of the 4.9 Mev. radiation. 


Coincidences (Arbitrary Units) 
ny 





3.0 3.5 4.0 45 5.0 5.5 
MEV. 


Fic. 3. Survey spectrum of carbon capture gamma-radiation. 


Fig. 4 shows the result of a detailed examination of the 4.9 Mev. radiation. 
The energy of this y-ray is 4.947 + 0.008 Mev. Correcting for the recoil of the 
carbon nucleus the corresponding transition energy is 4.948 + 0.008 Mev. and 
this must be the binding energy of a neutron to C™. Another value for this 
energy may be obtained by the addition of the following equations: 


C? + H? = N®4 n — 0.281 + 0.003 Mev. (2). 
N® = C84 Qmoc?+ 1.202 + 0.005 Mev. (14, 20). 
n + H! = H?+ 2.231 + 0.007 Mev. (1). 
n — H! = 0.782 + 0.002 Mev. (19). 


The result is: 
C2? + n» = C+ 4.954 + 0.010 Mev. 


Alternatively, the energy balance, Q, in the C” (d, p) C" reaction is 2.729 + 
0.009 Mev. (5). From this result and the binding energy of the deuteron we 


find: C2 +n = C® + 4,960 + 0.012 Mev. 


All three values are in good agreement. 
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Fic. 4. Detail of carbon radiation. The results of different experiments are indicated. 


Nitrogen 


The results of the survey measurements on the nitrogen capture radiations 
are contained in Fig. 5. The beryllium oxide spectrum is included in this figure 
for identification of the beryllium radiation. Aluminum radiation is present in 
all three spectra of Fig. 5 because all three compounds were contained in dural 
containers. The coincidence rates of the beryllium and aluminum peaks in the 
beryllium oxide spectrum are higher than those of the other two spectra 
because a thicker radiator was used. The difference between the peak counting 
rates of the nitrogen y-rays in the urea and beryllium nitride spectra can be 
accounted for quantitatively in terms of the difference in the nitrogen content 
of these materials. The nine radiations, lettered A to I in Fig. 5 are assigned 
to nitrogen. 


Figs. 6 to 10 show the results of more detailed studies of the nitrogen radia- 
tions obtained with beryllium nitride. Table I contains a summarv of the 
energies of these radiations and the relative intensities after correction of the 
peak coincidence counting rates for the energy sensitivity of the spectrometer.* 
The detail of the lower energy radiations is shown in Fig. 6. The peak at 4.70 
Mev. is due to the aluminum of the container. The feeble peak at I, however, 
is certainly due to nitrogen; its intensity relative to the 4.70 Mev. peak has 
been carefully measured both for urea and for beryllium nitride, the ratio of 
the two intensities being in the ratio of the numbers of nitrogen atoms con- 


* The relative intensities given in Table I are slightly different from those previously reported (15), 
minor corrections to the efficiency calculations having been made since. 
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Fic. 5. Survey spectrum of nitrogen capture gamma-radiation. The nitrogen radiations may 
be recognized by comparing the three spectra. 


tained by the two samples. The 4.95 Mev. y-ray due to capture by C” has 
been identified in a detailed study using urea; the coincidence counting rate 
was very low and only just enough to be distinguished above the statistical 
error of the measurements. The height of the peak was roughly equal to that 
calculated from the preceding carbon measurements. Fig. 6 shows some evi- 
dence for radiation at 5.11 Mev. on the tail of the peak H. A similar bump 
was observed on the tail of this peak in a detailed experiment using urea. In 
addition some indications have been found for radiations with energies of 5.66 
and 6.88 Mev. 


The y-ray F appears to be single and is shown in Fig. 7. The 7 Mev. radia- 
tions are shown in detail in Fig. 8; there is good evidence here for two radia- 
tions at 7.16 and 7.36 Mev. Since the neutron binding energy in N* is given 
by the energy of the y-ray A, viz. 10.82 Mev., the radiation at 7.36 Mev. might 
be expected to be in cascade with a radiation having an energy of 3.46 Mev. 
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Fic. 6. Detail of nitrogen radiations, G, H, I. 


TABLE I 
ENERGIES AND INTENSITIES OF NITROGEN ‘Y—-RAYS 


Intensity relative 


Y-ray Energy, Mev. oe that of A 
A 10.816 + 0.015 1.00 
B 9.156 + 0.030 0.09 
c 8.278 + 0.016 0.19 
D 7.356 + 0.012 0. 56 
E 7.164 + 0.010 0.19 
F 6.318 + 0.010 0.9 
G 5.554 + 0.010 1.5 
H 5.287 + 0.010 2.3 
I 4.485 + 0.010 0.8 
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Fic. 7. Detail of nitrogen radiation F. 
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If its intensity were equal to that of the 7.36 Mev. radiation, the peak coin- 
cidence counting rate should be 0.03 coincidence per minute after allowance 
is made for the energy sensitivity of the spectrometer. A search for this radia- 
tion was made between 3.0 and 3.75 Mev. but none was detected above a 
background of 0.05 coincidence per minute. 


r o 
7.1644 0.012 Mev. 7.3564 0.012 Mev. 


COINCIDENCES PER MIN. 
o fs a 


~~ 





69 70 71 72 73 74 MEV. 
Fic. 8. Detail of nitrogen radiations D and E. 


The 8.28 and 9.16 Mev. radiations are shown in Fig. 9. Both appear to be 
homogeneous. The ground state transition, i.e., the 10.82 Mev. radiation A, 
is shown in detail in Fig. 10. It is also homogeneous. 


The energies of the y-rays H and G, and those of I and F, add up to give 
the energy of the y-ray A, and therefore are in cascade. The binding energy 
of a neutron in N* may be obtained from the energies of the direct radiation 
A and the sum of the energies of these cascading radiations. Correcting for 
nuclear recoil and adding, these energies are 


From A 10.820 + 0.015 Mev. 

From H + G 10.843 + 0.014 Mev. 

From F + I 10.805 + 0.014 Mev. 
The average of these determinations is: 

N“ +n = N® + 10.823 + 0.012 Mev. 


The Q of the N" (d, p) N** reaction reported by Lauritsen (17) is 8.57 Mev. 
A more recent and probably more accurate determination is 8.61 + 0.02 
Mev. (18). The latter leads to a neutron binding energy in N" of 10.84 + 0.02 
Mev., which is in good agreement with the above value. 
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The neutron binding energies of C* and N*® may be used to calculate the 
values of certain mass differences which may be compared with mass spectro- 
graphic measurements. Taking into account the probable errors of the three 
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Fic. 9. Detail of nitrogen radiations B and C. 
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Fic. 10. Detail of nitrogen radiation A. 





results for the binding energy of C™, the mean binding energy is 4.953 + 0.006 ° 
Mev. From this we obtain: 
C#®H! — C8 = 44.80 + 0.06 x i0™4 mu. 
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Using our value for the binding energy of N"* we obtain: 
N“H! — N¥® = 107.85 + 0.11 X 10~ mu. 


The first of these two results is significantly higher than that obtained 
directly by Ewald, (10) viz., 44.10 + 0.08 x 10~* mu. Our measurement of 
the binding energy of a neutron in C” therefore confirms a discrepancy between 
the transmutation and mass spectrographie data indicated by the results of the 
measurement of the Q of the C” (d, p) C™ reaction. For the nitrogen difference 
Ewald obtained 107.76 + 0.20 X 107‘ mu, which is in agreement with our 


results. 
Discussion of the N® Levels 


The nitrogen radiations may be fitted to the level scheme shown in Fig. 11. 
With the exception of the new level at 9.16 Mev., this scheme is in agreement 


n'“ing)n® 


7.36 
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632 
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Fic. 11. Energy level diagram and capture gamma-ray decay scheme of N*. 


with the results of other authors. The 5.3 Mev. level was first observed by 
Cockcroft and Lewis (6) using the N“ (d, p) N™ reaction. Lauritsen (17) shows 
four levels below the binding energy of the neutron in N"; their energies are 
. §.29, 6.0, 7.2, and 8.2 Mev. The 9 Mev. level seems to have been observed as 
a proton group in the (d, p) reaction by Guggenheimer, Heitler, and Powell (12) 
with 6 Mev. deuterons. Most other observers, with deuterons of 1 Mev. or less, 
as might be expected, have not observed a proton group corresponding to this 
excitation. It is possible that the 9.16 Mev. y-ray does not arise from a tran- 
sition to the ground state from a level with this excitation energy; it can result 
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from a transition from the capturing state to a level at 1.66 Mev. A proton 
group corresponding to a level at 1.53 Mev. excitation has been observed by 
Holloway and Moore (13) using 1.07 Mev. deuterons. This group, however, 
has not been detected by other observers. 


Some of the nitrogen levels seem to possess a fine structure. Malm and 
Buechner (18) report the presence of two proton groups from the (d, p) re- 
action corresponding to excitations at 5.29 and 5.32 Mev. The resolving power 
of the present experiments has not been sufficient to separate such a close 
doublet. Wyly (21) has shown that both the 6.4 and the 7.3 Mev. levels are 
double. The energy of the deuterons in his experiments was insufficient to 
produce in its entirety the proton peak corresponding to the 8.2 Mev. level, 
but the work of Guggenheimer ef a/. shows that this level is probably double 
also. The present experiments, however, demonstrate the doublet character of 
the 7.3 Mev. level but give no evidence for doublet structure in the 6.4 and 8.2 
Mev. y-rays. The lack of such evidence is not surprising; one cannot be sure 
that a multiplet structure found by the use of charged particle reactions should 
appear in the y-ray spectrum. 


The y-rays produced by the N* (d, p) N® reaction have been studied by 
several observers with, unfortunately, little resolution. The first measurements 
to be made were those of Crane and his associates (8). y-Rays with energies 
of 8.2, 6.6, and 5.1 Mev. were observed by Crane, Halpern, and Oleson (9) by 
the study of Compton recoil electrons in a cloud chamber. These radiations 
clearly correspond to the excitation of the levels of Fig. 10. Gaerttner and 
Pardue (11), however, using a similar method, but studying the energies of 
electron pairs, found the 5.3 Mev. y-ray together with one at 7.2 Mev. Their 
results did not show the 6.6 or the 8.2 Mev. radiations mentioned above. This 
discrepancy has not been clarified. 


Intensities of the Nitrogen Radiations 


The intensities of the nitrogen radiations relative to that of the ground state 
transition, the y-ray A, are listed in Table I. In consideration of its high energy, 
the intensity of the y-ray A is surprisingly weak. Single particle nuclear models 
indicate that the parities of the ground states of N' and N" should be even 
and odd respectively. Since the spins of N" and N" are 1 and 1/2 respectively, 
the ground state transition must be electric dipole. However, A is no more 
intense than radiations G or I which have only about one half of its energy. 
It is worth noting that the intensity of the cascading radiations G and H are 
not equal. The y-ray H being more intense than G, it is clear that the 5.29 Mev. 
state is excited, in part, by transitions from higher levels; the radiations, 
corresponding to these transitions, however, are separately too weak, or possess 
energies too low to be observed. On the other hand, the y-rays F and I, which 
are also in cascade, are of nearly equal intensity. 
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Radiative Capture Cross Section for Nitrogen 


The total number of captures in nitrogen will be roughly proportional to the 
sum of the intensities of the radiations A, B, C, D, E, F, and H. The ratio of 
the cross sections for thermal neutron capture by nitrogen and by beryllium 
will be equal to the ratio of the sum of the nitrogen y-ray intensities (from 
beryllium nitride) to that of the beryllium y-ray, multiplied by 3/2, the ratio 
of the numbers of beryllium to nitrogen atoms in beryllium nitride. The cross 
section ratio is 18, whence it may be deduced, that if the cross section of 
beryllium is 9 millibarns, the radiative capture cross section of nitrogen is 160 
millibarns. If there is appreciable excitation of the 3.375 Mev. level in Be, the 
nitrogen cross section will be lower in proportion. The above value is therefore 
an upper limit for the cross section. 


This estimate of the nitrogen capture cross section leads to a rough value for 
the proton width of the capturing state. The total absorption cross section of 
nitrogen is 1.7 barns (7) and the penetrability of the Coulomb barrier for a 
proton of 600 kev. is about 1/27. If then the radiation width is of the order of 
a few volts (3), as would seem probable for light elements, the proton width 
without barrier is greater than this by a factor which is equal to the reciprocal 
of the penetrability multiplied by the ratio of the total absorption cross section 
to the radiative capture cross section. The proton width, therefore, is of the 
order of a kilovolt. This is a very low result. 
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THE ANGULAR DISTRIBUTION OF THE RADIATION AT THE 
INTERFACE OF TWO ADJOINING MEDIA' 


By S. CHANDRASEKHAR 


Abstract 


In this paper the results contained in a recent paper by B. Davison entitled 
“Angular Distribution of Neutrons at the Interface of Two Adjoining Media’”’ 
are derived by certain elementary methods based on certain simple invariances 
characterizing the problem. 


1. Ina recent paper published in this Journal, B. Davison (2) has considered 
a problem in the theory of radiative transfer which can be formulated as fol- 
lows: 


The half-space (r < 0) above a semi-infinite plane-parallel isotropically 
scattering atmosphere with an albedo a; for single scattering is filled with an 
isotropically scattering atmosphere, similarly stratified, and with an albedo w2 
different from w;. The equations of transfer valid for0 < 7 < ©» and0 >7r> 
— © are, respectively, 





rs 
palen) = T(x) -4e| T(r,n')du! (7 > 0) (1) 
Tv — 
+ 
and pO”) = T(r, 4) - bes | (r,u')du’ (7 < 0), (2) 


where 7 is the optical thickness measured from the interface of the two media 
inwards and uz is the cosine of the angle to the outward normal. 


It is known that Equation (1) admits an integral of the form 
I L eft 
sia i Mai carseat 3 
(su) = Le, (3) 
where Lp is a constant and k; is the positive root less than 1 of the character- 
istic equation 
w, log [(1 + ki)/(1 — ki)] = 2h: (4) 


Consider the solution of Equations (1) and (2) which has the asymptotic be- 
havior (3) for 7 — © and which vanishes for 7 — — o. The problem is: What 
is the angular distribution of the radiation at the interface, r = 0, of the two 
adjoining media? In the paper we have referred to, Davison has solved this 
problem by making use of the integral equations of the Schwarzschild—Milne 
type which govern the problem. However, we shall show that all of Davison’s 
results can be obtained by entirely elementary methods by appealing to certain 
invariances of the problem in the manner developed by the writer in recent 


1 Manuscript received August 31, 1950. 
Contribution from Yerkes Observatory, Williams Bay, Wisconsin, U.S.A. 
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years and summarized in his book (1) (hereafter referred to as Radiative 
Transfer). 


2. Let the law of diffuse reflection by a semi-infinite plane-parallel atmos- 
phere be expressed in terms of a scattering function S in the usual fashion 
(Radiative Transfer, p. 90). For an isotropically scattering atmosphere with 
an albedo oa the scattering function is given by (Radiative Transfer, pp. 97-98) 





S(u,u') = 0 *— Hw) Hv’), (5) 
btu 
where — uy’ refers to the incident direction and yu to the scattered direction and 
. 1 ? dy’ 6 
A(u) =1+ 45 Stu uw’) =F ) 
0 
satisfies the integral equation 
1 , A(u’) , 
H(u) = 1 + § on (x) —— dy’. (7) 
on+y 
An alternative form of Equation (7) which we shall find useful is 
1 * H(u') 
wi - ton] OD yy, (8) 
H(u) ee Jou + a 


3. Let I (r, + u) and IJ (7, — ») (0 < w <1) refer to the outward and the 
inward directed radiations at the level 7. Now the inward directed radiation 
I (r, — w) at a level tr > O must be reflected by the semi-infinite atmosphere 
below r according to the scattering function S,(y, u’), appropriate for an iso- 
tropically scattering atmosphere with an albedo w:. (We shall use a sub- 
script 1 to all functions referring to an atmosphere with an albedo w, and 
similarly a subscript 2 to all functions referring to an atmosphere with an 
albedo m2.) It isevident that this reflection of the radiation J (r, — u) by 
the atmosphere below 7 must account for the difference between I (7, + yu) 
and J (0, + u). Indeed, a little consideration shows that (cf. Radiative Trans- 
fer, p. 81, for the analogous formulation in the conservative case) 





—k,r fi 
I(r, + uw) = &" 10, +4) + < k Si(u, w’) I(t, — w)dw’ (r>0). (9) 
n 


The reason for the appearance of the factors e”’ and e~"” in Equation (9) is 
that we are considering a situation in which the asymptotic behavior at r = & is 


en 


ite 





I(r, uw) > Lo (10) 


Accordingly, there is an infinite source of radiation at great depths and the 
transformation tr — 7 + #, where ¢ is a constant, is equivalent to multiplying 
all outward intensities by e*“ and all the inward intensities by e 


—kit 
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Similarly, considering the inward directed radiation I(7, — yw) at a level 
7 <0, we must have 


—k,r fl 
I(r, — w) =e" 1(0, — u) hee | So(u, mu’) I(r, + udu’ (r< 0). (11) 
uw Jo 





We can obtain a pair of integral equations connecting 7(0, + uw) and 7(0, — yp) 
by differentiating Equations (9) and (11) and passing to the limits r = + 0 
and + = — 0, respectively. Thus, from Equation (9) we obtain 


; k , , ’ , 
pees # = kil(0, + ») — 3 | Si(u, w’)1(0, = 2’) du 
dr r=+0 2uJ0 


] 1 dI hi . , 
$s | Si(u,p') | Pe "| du’. 
2u 0 r 


dr =+0 


(12) 


On the other hand from the equation of transfer (1) (which is appropriate for 
t > 0) we conclude that 





[sie + 9) = 110, +») — 9:50) 
art r=+0O0 yp (13) 
iat jae * 2] tote tic tit: aue 
dr r=+0O 4p 


where we have written 


1 
J(0) = +] 10, a) du. (14) 


Substituting from Equations (13) in Equation (11) we obtain after some minor 
rearranging of the terms: 


(1 — kw) 1(0, + x) = wi J(0)| 1 + | 


1 : d , 
Sulu, 0) 4 | 
Be 


0 


1 du! (15) 
- +l Si(u, wu’) (1 + Rip’) I(0, — w’) ta . 
bu 
Similarly, from Equation (11) we obtain 
1 , 
(1+ bin) 10, — 2) = wJ(0)| 1 + +| Salus u’) = 
° # (16) 


1 du’ 
- | So(u, u’) (1 — kw’) 10, +n’) = - 
b 


Equations (15) and (16) are the mathematical expressions of the invariance of the 
angular distribution of the radiation at the interface adjoining the two media for 
arbitrary normal displacements of the interface. 


Now substituting for S; and S; according to Equations (5) and (6), we obtain 
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(1 — Rin) IO, + w) = w1 J(0)Ai(n) 


1 , 17 
— $e ntta(n) | AO) (1 + e100, = wr - 
on+y 


and 


(1 + ki w)I(O, — uw) = w2J(0)H2(u) 
18 
— Fouts) | ee (1 — ky’) TO, + w’)dy’. 7" 


Equations (17) and (18) represent a pair of integral equations governing 
I(0, + #) and I (0, — ht). 


Writing 
(1 — kyw)I(O, + w) = o1J(0) Fi (u) (19) 
and (1 + kiw)I(0, — w) = w2 (0) Fr(u), 
we observe that Equations (17) and (18) become 
1 ’ 
Fw) = Hw [1-4 omn| ae) Fau!)an' | 
Out 
t Fa(y") (20) 
and F(a) = Haw)|1— ein | 2 ryarae’ |. 
Opts 
From these equations it is at once evident that 
Hy(u) H2(p) 
F\(u) = —— and F,(u) = —— - (21) 
a id) “ Ayu) 
The angular distribution of the radiation at r = 0 is, therefore, given by 
10. + 9) © JO 
ix(u)(1 — ki) 
(22) 
Hyi(u) 


d I 0, = = weJ(0 —_—_——___———_ + 
” re Fd + kw) 


It can be verified that the foregoing solution is equivalent to that given by 
Davison (loc. cit. Equations [21] and [30)). 


4. For the solution (22) to be self-consistent it is necessary that 


1 
J(0) =4 I. (100, + ») + 10, — w)) du 
; Aiyi(u) [. H(z) \ 
= —_——_—_—__— d } ane 8 WA, a d . 
10){4e | ee dl a 


In other words, we must have 


ie oo Ay(u) bes| se dp = 1. (24) 
0 H2(u)(1 — = Hy(u)(1 + kw) 


(23) 
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The truth of this identity can be verified as follows: 


Substituting for 1/H2(u) according to Equation (8) in the first integral in 
Equation (24), we have 


1 [. ii (u) 


w es 
a oH —kw) 


: Ay(u) 1 1 i. Bae) ay . 
wi du ae — 7o2u Pere 


On the other hand, since —1/k; is a pole of Hi(u) (Radiative Transfer, pp. 123 
and 350), 


(25) 





i 


Law; | ty wt. (26) 
ol — ku 
The right-hand side of Equation (25) therefore becomes 
Hy (u)H2 (u’) ” 
l1—ioio c I. : aneesemsaserneestenntonmnene (an digs’, (27) 
oem (1 — kw) (u + v’) 


Now it can be readily verified* that 











1 ; pH s(n) TT 1 1 
en) Ne iy , (28) 
o(1 + bu)(u + wv’) 1 — by’ LAi(u’) =A, (1/8) 
where } is a constant. For b = —1/k:, this relation reduces to (cf. the remarks 
preceding Equation 26) 
1 
ben] ee ty (29) 
o(1 — kix)(u + wu’) (1 + kip’) Ai(u’) 


Accordingly the integration,over » in (27) can be performed and we are left 


with 
1 1 , 
al Ay(u) ae | Ai2(u’) du’ 30 

® lame ae “lies” fF 


This completes the verification of the identity (24). 


5. The net flux + F(0) of the radiation at r = 0 can also be directly eval- 
uated. Thus 


, , Ay(u)u 
2; (0, 201J(0) | ———————_. 
I. (0, + uw) udp = 2 wi J(0) o Hw 0) 


= 0 dp, —-1 1- D2 dy! 
ky s LI, r w{ war 0 ae 


1 
= 4 xof1 —{1-(1-o,)"”"} - pews | dutta) 








* By expressing u/[(1 + bu)(u + p’)] in partial fractions and using the equation satisfied by 
Hi, (y). 
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P Ai(u)u 1 Ayu) \ 
du ————_—_____ -— | d 
x1], "(ut u)0 — km) 0 rate 


_« = ote | ’ ' ENE EL ALS 


1 
et 1 =< 1/2 
Ai, (u’) ore ] 
- 4 x0)[ a — w:)¥? — (1 — w)¥? {1 — (1 — ws)*/*} 


1 AA(y’)p’ | 
+ $k | canner 
Fe Jeol + hw )Eie) 





(31) 


1 
- - (1 — wi)”9(1 — ws)27(0) + 2|/ 10, — »)udp. 


In the foregoing reductions we have made use of the various integral properties 
of the H-functions (Radiative Transfer, Chap. V.). We have thus proved that 


FO 4 — w;)/2(1 — ws)¥? 
70) i (1 — wi)¥?(1 — w2)¥?. (32) 


6. Finally, we may draw attention to a further invariance characterizing 
the problem. Considering the radiation at the interface r = 0, we observe that 
the inward directed radiation J (0, — yu) at r = 0 can be regarded as resulting 
from the reflection of the outward directed radiation J(0, + yu) by the semi- 
infinite atmosphere overlying r = 0. Thus: 


1 
mS | Sa(u, w’)TO, + udp’, 
~ (33) 
p'Hy(u’) 
pty’ 


On the other hand the reflection of (0, — 4) by the atmosphere below r = 0 
cannot account for all of J(0, + u); it can account for only that part of the 
outward directed direction at r = 0 which is in excess of what will be present 
in the absence of an atmosphere above r = 0. Now we know that when there 
is no atmosphere above r = 0 the angular distribution of the emergent radi- 
ation is given by (Radiative Transfer, p. 346) 


1 
or I (0, — ») = $02 H2 (u) I, I(0, + w’) dy’. 








Lo Aiy(u) 
I(0, = . 34 
Ot" Rap Pte on 
Accordingly 
Lo Ai(u) 1 I. ’ , , 
I(O, + p) — ————— ——- = — | Sile, p’)I(0, — d 
( #) Hil/k) 1— km 2udo i(u, w’)T( nu’) dp 


(35) 


1 \ 
= bos His) | Hie) 100, — »') de’. 
outh 
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It can be directly verified that with J(0, + yu) and J(0, — yu) given by Equations 
(22) the invariance (33) is identically satisfied. Moreover, the evaluation of the 
right-hand side of Equation (35) for J(0, — u) given by (22) leads to a deter- 
mination of the constant Ly in the asymptotic behavior of the solution at 
T= ©, We find 
Lo = o1J(0) Hi, (1/k:) ; (36) 
H2(1/k1) 
7. In conclusion it may be stated that there is no difficulty in extending the 
methods of this paper to more general laws of scattering (cf. Radiative Transfer, 
Chaps. VI, 1X, X, and Appendix III). 
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THE EFFECTS OF PLASTIC DEFORMATION ON MAGNETIC 
PROPERTIES OF POLYCRYSTALLINE METALS! 


By Ursuta M. Martius 


Abstract 


Most of the laws of the ferromagnetic behavior of metals and alloys were 
derived from studies of single crystals. In order to apply these laws to poly- 
crystalline metals the effect of the crystallographic composition of the poly- 
crystalline body has to be taken into account. For this purpose the concept of a 
“closed flux shell’ as a fundamental quality of polycrystalline metal is put 
forward on the basis of the following trend of thought: 

One of the most important factors determining the size and arrangement of 
ferromagnetic domains in a single crystal is the tendency to close the magnetic 
flux throughout the specimen by a suitable arrangement of these domains. The 
size of the crystal defines the limits within which the flux can be closed. Ina 
polycrystalline metal, a “‘crystal’’ is defined as a region of one crystallographic 
orientation separated by crystal boundaries from neighboring regions of differ- 
ent orientation. The corresponding unit, as far as ferromagnetic properties 
are concerned, will be a closed flux shell which may contain many domains but 
is always a region of closed magnetic flux. Neighboring regions, having differ- 
ent crystallographic orientation, will each also consist of a closed magnetic flux 
shell. Crystal boundaries are regarded as lattice discontinuities which separate 
regions of different orientation and at the same time separate different closed 
flux shells. 

The consequences of this model and experimental evidence supporting it are 
discussed. An explanation of the effect of plastic deformation on magnetic 
properties is suggested by considering the changes in the size of the closed flux 
shell which will occur during plastic deformation. 


Recently, both the theory of mechanical behavior and the theory of ferro- 
magnetic properties of metals have been advanced considerably. At present 
it seems worth while to investigate the effect of mechanical changes in metals 
on their magnetic properties with the object of using the observed changes in 
magnetic properties to obtain a more thorough understanding of the mechanism 
of deformation and failure of metals. 


The Ontario Research Foundation has been engaged for some time in such 
an investigation. In utilizing present ferromagnetic theory to explain the 
observed behavior of magnetic materials during deformation, it became quite 
apparent that the magnetic properties were not changing in our polycrystalline 
samples in any way that could be explained fully by the domain theory as we 
know it. This theory has been derived mainly from the study of single 
crystals. .It became evident that the introduction of a new unit of magnetic 
structure was necessary in our consideration of the behavior of polycrystalline 
metals. 


- The conception of a ‘‘closed flux shell" as the unit of magnetic structure in 
a polycrystalline metal is put forward in this paper. This conception has 


1 Manuscript received July 11, 1950. 
Contribution from the Department of Engineering and Metallurgy, Ontario Research 
Foundation, Toronto, Ont. Presented at the fifth Annual Congress of the ‘anadian Association 
of Physicists, Mc Master University, Hamilton, Ont., May, 1950. 
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proved useful in visualizing and predicting the magnetic behavior of poly- 
crystalline metals under the influence of external stress, particularly in the 
range where plastic deformation of the metal occurs. In order to discuss the 
consequences of this conception, it will be necessary to give a brief account of 
the present status of the dislocation theory of plastic deformation and of the 
domain theory of ferromagnetism. 


The Dislocation Theory and Plastic Deformation 


The current ideas on the strength of solids are based mainly on the dis- 
location theory (5). 


According to this theory, every real crystal contains a certain number of 
dislocations. A dislocation is a linear lattice defect which exists in a crystal 
when two atomic planes are displaced in respect to each other over part of 
their surface. There are two main types of dislocations, namely ‘‘edge” and 
‘‘screw’’ dislocations, the difference depending on the direction of displacement. 


The dislocation is characterized by its mobility along the glide plane 
and the fact that such movement causes the one half-crystal on either side of 
the glide plane to be displaced in the slip direction. Edge and screw dislo- 
cations differ only in the direction of the dislocation axis relative to the slip 
direction. These specific properties distinguish dislocations sharply from 
other types of lattice defects. A dislocation as a lattice discontinuity cannot 
end within the crystal but continues on to the crystal boundaries or else forms 
a closed loop. It is possible to treat grain boundaries as arrays of disloca- 
tions, as Shockley and Read have done recently (10). 


Because of the atomic displacements, which define them, dislocations are 
centers of internal stress. The stress falls off inversely with the distance from 
the center of the dislocation. Through this stress field a dislocation can 
interact with other dislocations, with crystal boundaries, foreign atoms and 
irregularities of the lattice. In many cases the mechanism of interaction can 
be predicted. If, for instance, an atom in an interstitial position causes a 
local lattice expansion, it will migrate under the action of the stress field of a 
dislocation to the locally expanded region at the crystal around the dislocation. 
Atoms causing compression will consequently concentrate in the compres- 
sion zone. 


It is the basic conception of the theory that plastic deformation in crystalline 
solids takes place through the motion of dislocations. This idea has proved 
extremely successful because, in connection with the established properties of 
dislocations, it affords an explanation of creep and slip, of the existence of an 
elastic limit and its dependence on preworking and other properties of the 
metal. 


If an external force is sufficiently strong, it will cause a dislocation to move. 
That is equivalent to saying that the forces will cause slipping of one part of 
a crystal with respect to another. Considering the details of this mechanism, 
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it is postulated that new dislocations can be created in the same slip plane 
which then move under the same external forces. This process cannot be 
extended indefinitely, since resistance to further deformation increases with the 
amount of plastic deformation a metal has already undergone. This can only 
mean that it becomes more difficult to produce or move a dislocation. The 
forces necessary to move a dislocation increase with increasing numbers of 
dislocations. The stress fields of the dislocations interact and this interaction 
limits the dislocation density. Eventually the positions of the dislocations 
with respect to each other will be fixed in some regular pattern. On this basis 
it is possible to understand why the deformed metal has a higher elastic limit 
and why the resistance to further deformation increases with the amount of 
previous deformation. These new conceptions will have considerable bearing 
on our interpretation of the ferromagnetic properties of metals. 


Domain Theory of Ferromagnetism 


The main achievement in the recent work on problems of ferromagnetism 
has been the explanation and prediction of domain structures under given 
experimental conditions (7). 


It has long been realized that a ferromagnetic body is composed of “‘domains” 
magnetized to saturation and separated by transition layers called Bloch walls. 
The domains are so oriented that in the absence of an external field no over-all 
magnetization occurs. The observed magnetization process can always be 
explained by changes in domain size and changes in direction of magnetization 
of a domain with respect to the field. 


The thermodynamical reason for the formation of domains is that the free 
energy of the specimen can thus be reduced considerably. On the basis of 
the conception of minimum energy it is possible to derive domain structures 
for a given crystal. These structures have been verified experimentally by 
Williams and coworkers who also succeeded in showing the change in domain 
size caused by the action of an external field. These experiments are the most 
important confirmation of the basic ideas of the domain theory (11, 12). 


It is known today that there are two characteristic kinds of domains: interior 
domains and surface domains. The domains in the interior are fairly large 
and of simple structure. Their directions of magnetization are antiparallel (9). 
The surface domains are more complicated and their structure is mainly de- 
termined by the surface energy of the specimen. Their arrangement is always 
such that they close the magnetic flux throughout the specimen in order to 
avoid ‘‘free magnetic charges” on the surface of the sample. 


The idea that the tendency to establish a closed flux structure is one of the 
most important factors determining the domain structure was first put forward 
by Landau and Lifshitz, then taken up and extended by Néel and others (9). 


The actual size and arrangement of the domains depends largely on size, 
shape, and condition of the sample and the crystallographic properties of the 
metal. 
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Ferromagnetic single crystals have strong directional properties. They 
show certain preferred crystallographic directions along which they can be 
magnetized with a minimum amount of external energy. Considerable addi- 
tional energy is required to turn the magnetization from this ‘‘easy"’ direction 
into a ‘‘hard’’ one. This anisotropy energy relates the ferromagnetic to the 
crystallographic properties. It is also the root of the phenomenon of ‘‘magneto- 
striction’. Positive magnetostriction means that during magnetization the 
linear dimensions of a crystal increase in the direction of the applied field and 
that, conversely, under stress the magnetization increases in the direction of 
the applied tension. By virtue of this fundamental relationship, it is possible 
to influence the magnetization process by external stresses. 


There is, however, also an indirect effect of the stress. Every external stress 
will cause strains in the lattice and these deviations from perfect regularity will 
influence the domain behavior. Roughly speaking, the internal strains increase 
the resistance towards movement of the Bloch walls and cause changes in the 
anisotropy conditions. 


These considerations are valid only for stresses below the elastic limit of the 
metal. 


Magnetic Structure of Polycrystalline Metals 


The laws of the behavior of ferromagnetic domains outlined above are 
derived from investigations of ‘single crystals. For a polycrystalline metal 
without special texture, the same considerations should be valid. However, 
quite generally, the properties of a polycrystalline metal are not just the sum 
of the corresponding properties of all single crystals, suitably averaged. There 
are essential differences: for instance, the influence of the crystal surfaces and 
of the presence of the neighboring crystals will enter the picture and these 
effects will become increasingly prominent with decreasing grain size. 


A polycrystalline metal can be described as an assembly of small crystals 
oriented at random and separated by crystal boundaries. A crystal would 
then be defined as a region of one crystallographic orientation separated from 
the neighboring similar regions with different orientation by a crystal boundary. 
If this polycrystalline metal is ferromagnetic, the size and shape of the crystals 
will strongly influence the domain structure. As outlined for the case of 
single crystals, the presence oi the crystal surface and the fact of the finite size 
of the crystal are the reasons for the formation of the surface domains which 
tend to close the magnetic flux throughout the specimen. Thus, as far as the 
ferromagnetic properties are concerned we regard a “crystal” as a region 
which may contain many domains but has a single closed magnetic flux 


structure. 
The neighboring ‘‘crystal”, again a region of closed magnetic flux in itself, 


will have a different crystallographic orientation, and the whole polycrystalline 
metal can be pictured as an arrangement of these closed flux regions. 
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It is a purpose of this paper to introduce this concept of ‘‘closed flux shells’’ 
as a fundamental quality of polycrystalline ferromagnetic metals. As long 
as a region in the polycrystalline metal has the same crystallographic orien- 
tation, it will also form one closed flux shell. Any lattice discontinuity that 
acts as a boundary in the crystallographic sense by separating regions of 
different orientation will also act as a surface in the ferromagnetic picture by 
giving rise to the formation of surface domains and dividing closed flux regions. 


The polycrystalline metal would then be composed of closed flux shells 
separated by shell surfaces. Each shell would consist of an arrangement of 
domains. The domains within the shells are separated by Bloch walls. 


There is a fundamental difference between Bloch walls and shell surfaces. 
The Bloch walls are transition layers between regions of different magnet- 
ization, that is, different spin orientation (but it is assumed that the angle 
between neighboring spins is small and that there is no change in lattice 
parameter). The shell surfaces, on the other hand, separate regions of different 
crystallographic orientations and, at the same time, regions belonging to 
different magnetic ‘‘circuits’. Bloch walls can be displaced by an external 
magnetic field, whereas the field cannot act on the shell surfaces as such. Its 
only effect on a surface is an indirect one, caused by the change in the surface 
conditions through the displacement of the Bloch walls. 


The behavior of polycrystalline ferromagnetic metals under external forces 
will then be the superposition of two factors; (a) the actiou of the external 
forces on the single closed flux shell (where the laws of the behavior of single 
crystals are completely applicable), and (6) the action of the external forces 
on the arrangement of the shells. 


In the case of an external magnetic field, this may seem trivial. In this case 
it is always possible to take the action on the arrangement into account by 
considering the change in the internal demagnetization field, but even here the 
closed flux shell model points to some further conclusions: 


Within a closed flux shell the size and the distribution of the domains de- 
pends on the shape and size of the shell and on the conditions of the shell 
surface for any given metal. With decreasing shell size, which corresponds 
to decreasing grain size, the surface domains will become more and more 
predominant. They are separated by 90° walls so that with decreasing shell 
size the influence of the 90° walls will increase. This should affect the values 
of the coercive force and the hysteresis losses. 


Effect of External Stress on Magnetic Properties 


If external stress is acting on the polycrystalline metal, the influence of the 
two different factors of the magnetic structure becomes even more obvious. 
As long as the stress is, for every crystal, below the elastic limit, there will be 
only the action of the stress on the domains within the shell in the way outlined 
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above for single crystals. Experimental evidence for this was again provided 
by powder pattern technique (11, 12). However, as soon as the stress exceeds 
the elastic limit in any grain, a new mechanism starts because now the mechan- 
ical state of the sample changes. Up to this point we could regard the number 
and size of closed flux shells in a given specimen as constant. These will 
depend on pretreatment and chemical composition of the metal, but will not 
change during the various magnetization processes. As soon as plastic de- 
formation starts, this assumption is valid no longer. 


A crystal, after deformation, has undergone displacement of the part above 
the slip plane with respect to the part below. During slip new dislocations can 
be created in the slip plane and after slip the parts of the crystal now separated 
by an array of dislocations have a tilt of approximately 10 to 15 min. of arc. 
These changes in the mechanical state also change the magnetic properties. 


There are numerous experimental results showing how the ferromagnetic 
behavior of a metal under stress changes abruptly when the stress exceeds 
the elastic limit: This is much easier to understand with our model in mind. 


The direction of slip is for most metals a direction of hard magnetization. The 
dislocations in the slip plane, the aggregation of foreign atoms around these 
dislocations, the angular displacement of one crystal part with respect to its 
neighbor, all these will make the slip plane and its vicinity not only a new 
boundary in the crystallographic sense but also a new shell surface in the 
ferromagnetic sense. The original closed flux shell will split into two new 
shells and the slip plane will be the new shell surface. 


The basic assumption of this model is that the lattice discontinuities, which 
are formed in the course of plastic deformation, not only separate two regions 
of different crystallographic orientation, but also two different closed flux 
regions. 


In support of this assumption, there is first of all the experimental evidence 
afforded by numerous observations of the change in magnetic properties of 
metals under stress. 


There are two principle types of experiments being carried out to relate the 
changes in magnetic properties of a metal to its mechanical characteristics: 


The experiments of the first type measure the magnetization curve of the 
sample at a given constant stress. The observed change in magnetic proper- 
ties is then related to the changes in the mechanical state. The results of 
these types of experiments are obtained by measuring magnetic properties 
under a known load that remains constant during the individual measure- 
ment. That must be kept in mind when studying these data, even where 
they are eventually plotted as magnetic properties vs. stress. The experiments 
of the second type are of another nature. Here an attempt is made to record 
structure sensitive properties during the course of loading in order to follow 
the process of mechanical changes. 
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The structure sensitive property recorded in Type II is usually the effective. 
coercive force obtained by measuring the hysteresis loop at low frequencies 
and low fields, by recording the core losses in a sample in a very weak field 
and at frequencies of several kilocycles or by similar setups. 


Type I contains the classical experiments on the influence of stress on 
magnetic properties by Ewing, Becker, and others (1). 


The conditions of these experiments make it difficult to follow the changes 
near or above the elastic limit of the metal, but, even so, they indicate the 
complete change in magnetic behavior of a plastically deformed metal. Fig. 1 
shows the change in induction under stress and illustrates the rapid change in 
the neighborhood of the elastic limit (6). 
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Fic. 1. Change in magnetic induction of steel under tension. M. F. Fischer, 1928 (6). 


Figs. 2 and 3 show the results of experiments somewhat in between Types I 
and II (8). Here the change in induction is recorded for a specimen under 
repeated constant loading. The induction remains constant as long as the 
load is below the elastic limit. The slight initial increase is due to relief of 
internal stress during the first period of loading. As soon as the stress reaches 
the elastic limit of the metal, the induction increases abruptly. If repeated 
loading in this range is applied for a short enough interval so that failure does 
not occur, the magnetic characteristics of the metals change completely after 
removal of the load, as shown in Fig, 3. 


- Fig. 4 gives a record of an experiment representative for Type II (4). In 
this case the diameter of the minor hysteresis loop under low frequencies 
(125 c.p.s.) and low fields was measured for steel samples. The samples were 
loaded above the elastic limit (5). The load was released at the same speed, 
and the loading process was then repeated in the same way, The solid curve 
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corresponds to loading, the dotted curve to. unloading. Under the given ex- 
perimental conditions, the X-axis intercept of the hysteresis loop is inversely 
proportional to the effective coercive force. This means that the effective 
coercive force of the original sample remains constant till the load reaches 
about 5000 lb. Here too, initial increases at low loads occur sometimes, owing 
to the relief of interal stresses. 
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Fic. 2. Change in magnetic induction during repeated loading. E. Moench, 1940 (8). 


The data obtained at the first low loading depend on the individual mechan- 
ical state of the sample, whereas the remainder of the curve is rather inde- 
pendent of individual variations and representative for the general behavior of 
steel under stress. A further increase in load gives a sharp increase in effective 
coercive force at a stress level which is below the elastic limit. This increase 
is due to mechanical changes initiated by the start of plastic deformation in 
the individual crystals. Relief of the load brings first a slight decrease in the 
effective coercive force because of the mechanical recovery of the metal and 
the corresponding decrease of the internal demagnetization field. This process 
is completed around 5000 Ib. and the effective coercive force, now larger than 
in the undeformed metal, remains constant till the load has decreased to about 
3000 Ib. Further removal of the load brings a sharp increase in effective 
coercive force. 


This increase, which can be overcome by renewed loading, can be explained 
as follows: In our model, plastic deformation means that some of the original 
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closed flux shells are split up into new smaller shells. The creation of new 
shell surfaces results in the formation of new surface domains which have 90° 
walls toward the basic domains. These basic domains are separated by 180° 
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Fic. 3. Change in magnetic induction during repeated loading. Load briefly increased above 
the elastic limit. E. Moench, 1940 (8). 
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Fic. 4. Width of minor hysteresis loop under increasing tension (8). 


walls. As far as the action of the stress is concerned, magnetic domains sepa- 
rated by 180° walls are symmetrical toward the direction of stress and their 
volume will not change under external stresses. Domains separated by 90° 
walls, however, are affected by external stresses. Those favorably oriented 
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toward the direction of stress will grow by wall movement at the expense of 
those unfavorably oriented. 


The effective coercive force which we measure by following the hysteresis 
loop can be pictured as being proportional to the resistance against wall dis- 
placement. This resistance is decreased by the action of the stress on the 
90° walls. When all 90° domains are oriented as favorably to the direction 
of stress as it is possible to achieve by wall movement, the additional effect of 
the stress does not enter the picture any more. This stage is reached at a 
load of about 3000 Ib. and from there on Curve 2 follows Curve 1. Repeated 
or decreased loading does not change this pattern any further. 


RELATIVE CORE LOSSES 





STRESS, P.S.1. 


Fic. 5. Relative core losses under increasing tension (weak field, 5 kc.). P. E. Cavanagh, 


1946 (4)e 


Fig. 5 gives similar direct evidence of the changes in magnetic properties 
under plastic deformation (3). They show recordings of magnetic core losses 
under stress obtained in very weak fields of a frequency at 7 kc. (Curve A). 
The losses increase under stress as expected and then start to decrease below 
the elastic limit. This decrease in core losses corresponds to the decrease 
in effective coercive force during plastic deformation in Fig. 4. After plastic 
deformation the over-all core losses have increased (Curve B). 


The effects of lattice discontinuities on the magnetic properties will be more 
pronounced the larger the anisotropy energy of the metal. This is confirmed 
by an observation of Bozorth and coworkers who reported domains in a 
material of low anisotropy energy going right across grain boundaries as an 
exceptional case (2).* Williams showed the appearance of new lines in the 
powder pattern of plastically deformed permalloy crystals. These lines re- 
mained even after repeated magnetization and demagnetization and reappli- 
cation of the colloid in the same place on the crystal. They could not be 
observed without the application of the colloid. Williams relates these lines 


* The Ontario Research Foundation is indebted to Dr. Bozorth for this private communication. 
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to the plastic deformation of the crystal and suggests the use of the powder 
technique to observe early stages of plastic deformation (11, 12). 


It seems rather obvious that these lines are new shell surfaces created by 
plastic deformation at the slip plane, since they cannot be displaced by a 
magnetic field. Shell surfaces separate different regions of closed flux, so it 
should be possible to detect them by means of the powder technique. 


Néel has derived the structure of domains around impurities or cavities in 
ferromagnetic metal, and Williams has obtained powder pictures which verified 
Néel’s prediction. There seems to be little doubt that, if these impurities or 
cavities or their arrangements have linear dimensions comparable to the 
dimensions of the domains, this would lead to the formation of a new shell 
surface, surface domains and two separate closed flux configurations. We are 
going to carry out more quantitive investigations in this line. 


The concept of a closed flux shell outlined above seems to be consistent with 
the present experimental material and the theoretical picture of the ferro- 
magnetic as well as of the mechanical properties of the metal. It is suggested 
as an explanation of the discontinuities of the ferromagnetic behavior of 
plastically deformed metals. It helps us to visualize the effects of plastic 
deformation on magnetic properties and attempts to provide means for con- 
sistent application of the results of investigations on single crystals to poly- 
crystalline ferromagnetic metals. 


The necessary conditions of crystallographic orientation and anisotropy 
energy which will cause a given slip plane to act as a shell boundary are going 
to be derived quantitatively in a forthcoming paper. 
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THE ANGULAR CORRELATION OF ANNIHILATION RADIATION! 


By P. E. ARGYLE AND J. B. WARREN 


Abstract 


The angular distribution of annihilation radiation from copper 64 has been 
measured with a pair of scintillation counters in coincidence. Coincidences were 
found when the counters and source were not quite collinear. The relation 
between the noncollinear coincidence rate and angle of deviation A could be 
closely approximated to by an exponential of the form e—/», where Ay = (3.60 + 
0.04) 10-* radian. Interpreting this departure from antiparallelism of the two 
photons as arising from the residual motion of the annihilating pair, the average 
momentum of their mass center in copper is found to be (7. 20 mc.) X 10-3 which 
corresponds to an electron energy of 13.2 ev., if the positron is assumed to be 
thermalized before annihilation. 


Introduction 


By 1934 the general predictions of the Dirac theory as to the properties of 
positive electrons had been verified. Since then only a few of the finer details 
of the behavior of positrons have been investigated. One of these, the exact- 
ness of the collinearity of the two annihilation quanta, was investigated by 
Beringer and Montgomery in 1942 (1) who found no evidence for any lack of 
collinearity within the angular resolution of their experiment. It seemed 
worth while to repeat this experiment using scintillation counters and a 
stronger source, to ascertain whether all positrons are thermalized before 
annihilation and to obtain the order of magnitude of the momentum of the 
center of gravity of the pair at the moment of annihilation. 


Our measurements are similar to those recently described by DeBenedetti 
et al. (2) and our results are in agreement with theirs. 


Experimental Arrangement 


Two scintillation counters, one fixed and one mounted on a traveling micro- 
scope carriage, movable in a horizontal direction were placed on either side of, 
and 200 cm. from, a strong copper 64 source. The source, about 40 millicuries in 
strength at the start of the experiment, was a piece of copper 6 by 3 by 0.6 mm. 
irradiated in the Chalk River pile. The length of the source parallel to the 
line joining the counters was limited to 6 mm. to avoid attenuation of the 
annihilation radiation but this was sufficient to stop most of the positrons by 
self-absorption. No heavy absorbing material was placed over the copper to 
stop positrons emerging, since the thickness required would be such that more 
annihilation radiation would be absorbed than produced and the angular 
resolution of the apparatus would be reduced owing to the increase in source 
width. The copper was supported in a thin aluminum frame. The whole 


1 Manuscript received August 28, 1950. 
Contribution from the Department of Physics, University of British Columbia, Van- 


couver, B.C. 
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source was small in comparison with the detector size so that corrections for 
source size were unnecessary. 


Two clear rectangular blocks of naphthalene 2 by 1 by 1 in. were used as 
phosphors with a pair of RCA 5819 photomultipliers. The tubes were run at 
a relatively low potential (900 to 1000 volts), obtained from a battery refer- 
enced stabilized supply, and were cooled via their bases with solid carbon 
dioxide. This simple cooling reduced the noise background count, with the 
discriminator setting used, to about 40/sec., a negligible amount. 


The output pulses were fed via head and main amplifiers (rise time of 
0.2 ywsec.) into a conventional mixer of resolving time 0.30 usec. and with 
recovery time set at 20 usec., the longest dead time in the electronics. Care 
was taken to minimize overshooting in the amplifiers by the use of diode pumps. 
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Fic..1. Plot of coincidence rate against displacement. 
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Fic. 2. Logarithmic plot of tail of coincidence rate-displacement curve. 
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Procedure 


After careful alignment of the detectors and source the coincidence rate and 
the rates in each channel were recorded as a function of the horizontal displace- 
ment of the movable detector. Corrections for accidentals, dead time losses, 
and source decay were made to the observed coincidence rate in that order. 
Fig. 1 gives the general shape of the curve obtained and Fig. 2 a logarithmic 
plot of the exponential tail. The method of least squares was used to obtain 
count y 


the best slope, each point being weighted by the factor (__<omt _ 
™ ™ wine . standard deviation 


Results 


It was found that the tail of the distribution curve could be accurately 
represented by an exponential form, e~”™*, where Ao = (3.60 + 0.04) 1073 
radian and ) is the angle of departure from collinearity of the two photons. 
This value includes the effect of the finite angles subtended in the two planes 
by the crystal at the source. The reduction of this measured curve to a 
momentum distribution of the center of mass of the annihilating pairs was 
done by the method of DeBenedetti et al. (2). The mean center of mass 
momentum for annihilating pairs in copper corresponding to this value of Xo is 


Peu = 0.98 X oo which is equivalent to an electron energy of 13 ev. 


No corrections were applied in this reduction of the measured curve for the 
effect of small angle Compton scattering in the source nor for the effect of the 
finite detector sizes on the shape of the distribution. The Compton scattering 
would raise the counting rate at the outer end of the tail by about 1% which 
would not significantly change the slope of the curve. The effect of the size 
of the naphthalene blocks is to flatten the distribution curve, which might 
decrease the slope of the tail by perhaps as much as 5%. 


Discussion 
The value obtained lies within the probable experimental error of Dumond (4) 


who gives a value of 1.1 ns for p., deduced from the wave length profile of 
of 


the annihilation radiation. It compares with DeBenedeiti’s value of 1.2 for 
gold. 


The result is in agreement with the view that the majority of the positrons 
are thermalized and then annihilate with the outer and conduction electrons 
in the copper lattice. Measurements of this type in which positrons are slowed 
down in different materials would thus yield data as to the energies of the outer 
electrons, information of some interest in solid state theory. Thus far this 
technique does not indicate whether the theoretical predictions (5) regarding 
the annihilation of a significant fraction of positrons in rapid motion are 
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correct (3), since the maximum center of mass momentum observed corres- 
ponded to an energy of about 500 ev. 
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THE SCATTERING LENGTHS OF THE DEUTERON! 
By D. G. Hurst anp N. Z. ALcocK 


Abstract 


The angular variation of the scattering of neutrons of 0.0724 ev. energy by 
deuterium gas has been measured. Comparison with theoretical calculations of 
this variation gives a ratio of scattering lengths: 


a3 /a3 = 0.12 + 0.04 or 3.2 + 0.3. 


In combination with the known total scattering cross section of the free deuteron 


of 3.44 + 0.06 barns this gives for the scattering lengths the two sets of values: 
= 0.07 0. 826 0.012 
oe en \ or ? P \ xX 10-* cm. 
az = 0.638 + 0.006f 0.26 + 0.02 f 


I. Introduction 


Information about molecular structures and nuclear properties may be 
obtained from measurements of slow neutron scattering by gases. Preliminary 
measurements at Chalk River on oxygen and carbon dioxide gases (1) showed 
the feasibility of measuring the scattered neutron intensity as a function of 
angle. Consequently further work was undertaken with a new spectrometer; 
and in the present paper, a precise measurement of the angular variation of 
scattering by deuterium is reported. 


The scattering and absorption of slow neutrons by free deuterons is com- 
pletely defined by two scattering lengths (5), a3, and a; corresponding to the 
two possible spins 3/2 and 1/2 of the compound nucleus. In general, scattering 
lengths are complex, but when absorption is negligible, as is the case with 
deuterium, they are real. Accurate values of these scattering lengths are of 
interest for they are quantities with which theories of nuclear forces may be 
compared. 


The structure of the deuterium molecule is well known. The expected scatter- 
ing has been calculated by the methods of Hamermesh and Schwinger (8, 9) 
and Spiers (16, 17), and in addition corrected for the ground state vibration. 
Comparison with the measured scattering gives the value of (a3 — a})°(2a3 + 
a,)~*. This is a quadratic function of the ratio aj/a3 and has two roots. The 
scattering lengths corresponding to each of these roots are obtained by com- 
bining each value of a;/a3 with o7, the measured total cross section of free 
deuterons (6) which is equal to 4x (3a3? + 4a)”). This leads to four values 
of each scattering length. Two of these are eliminated by the crystal diffraction 
result, discussed in the next paragraph, that 2a3+ ay is positive. The remaining 
ambiguity could be resolved only by experiments with polarised neutrons (8). 


Two previous investigations in principle give the scattering lengths, but 
owing to lack of precision the authors did not carry the analysis to this end. 
1 Manuscript received August 22, 1950. 


Contribution from the Atomic Energy Project, National Research Council of Canada, 
Chalk River, Ontario, Canada. Issued as N.R.C. No. 2284. 
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Shull e¢ al. (15) using the diffracted intensity from powdered sodium deuteride, 
and later thorium deuteride (Wollan e¢ al. (18)) obtained the coherent scattering 
length fp = $ (a3 + 404) and showed that it is positive. A measured value 
of fp combined with a measured o7 gives the scattering lengths. Fermi and 
Marshall (6) measured the transmission cross section of deuterium gas for very 
low energy neutrons, and the total cross section of the free deuteron. From 
the results, they gave rough limits for the ratio a}/az and showed that it is 
positive. In Section V1I(b) below, a comparison is made of these methods with 
the present work. 
II. Theory 

(a) The Deuterium Molecule 

The deuterium molecule is diatomic, with identical nuclei having spin 1. 
Owing to the symmetry properties of the molecule, there are two parities, 
depending on whether J, the quantum number of the rotational state, is even 
or odd. States of even J are called orthodeuterium, and odd J paradeuterium. 


The rotational energy of the state of quantum number J is 


(h/2x)? J(J +1) 
E; omn? +) Ey = 
where h = Planck’s constant, 
M = mass of the hydrogen atom, 
r. = separation of the nuclei in a rigid model of the molecule. 


At a temperature T the Jth state has statistical weight 
W, G(2J + 1) exp (—E,/k,T) 


G = 6 for J even, 
= 3 for J odd, 
k, = Boltzmann’s constant. 


The gas contains 66.7% orthodeuterium at room temperature, 70.3% at liquid 
nitrogen temperature. Transitions between the two parities are very rare and 
are of importance only under the action of catalysts. When the gas is cooled 
from room temperature the new equilibrium is set up slowly. 


In addition to rotational states, the molecule possesses vibrational states. 
The excitation energy of the vibrational states is so large that in the scattering 
gas only the ground state occurs. The wave function of the ground state extends 
over about one tenth of the internuclear separation. This modifies the neutron 
scattering in that the scattering material is not confined to a thin spherical 
shell as on the “rigid’’ model. 


The Maxwellian distribution of the velocities of the molecules influences the 
scattering by affecting the relative momentum of the neutron and molecule. 
At room temperature, the motion of molecules of average velocity causes a 
difference in the effective neutron wave length of about + 30% for the ex- 
tremes of motion away from and towards the neutron. At lower temperatures, 
the spread of effective wave lengths is correspondingly reduced. For this 
reason, and to have few rotational states excited, the experiments were per- 
formed at low temperature. , 
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(b) Outline of the Calculations 


The basic formulas for the scattering of neutrons by deuterium were derived 
by Hamermesh and Schwinger (8) who were concerned mainly with the be- 
havior of the total cross section. To put their results for the differential cross 
sections into suitable form to compare with our measurements required that 
the cross sections be transformed from the center of mass system to the labor- 
atory system, and averaged over the Maxwell distribution of molecular 
velocities. 


Spiers has done this (16) and gives the cross sections for molecules at rest 
for scattering angles between 0° and 180° and the averaged cross sections for 
molecules in a gas at a temperature of 90°K. (liquid oxygen temperature) at 
scattering angles of 30°, 60°, 90°, 120°. From his report the region of scattering 
angles between 10° and 70° was chosen as the most promising for an experi- 
mental determination of the ratio of scattering lengths. Accordingly a second 
set of cross sections was computed for a gas temperature of 77.4°K. at angles 
of 10°, 30°, 50°, 70°. The temperature 77.4°K. is the temperature of liquid 
nitrogen, the refrigerant used in the experiment. 


There is a partial cross section for each transition J — J’. These partial 
cross sections have to be summed over all J’ and J, the statistical weights being 
applied in the sum over J. The result is a differential cross section o(@) expressed 
as a linear combination of two quantities ¢4 and og which can be calculated 
without a knowledge of the scattering lengths. Thus 


o(¢) = Aca($) + Boz(9), (1) 


¢ being the angle of scattering in the laboratory system. The scattering lengths 
occur only in the coefficients A and B. 


The measured scattering s(¢) should be proportional to o(¢) and we may 


write c*(¢) = Ao,+ Bog which transforms to s() = (Ao,/og + B)/c, where 


oB 
c is a constant of the apparatus. A plot of s(¢)/og against o4/o, should be a 
straight line cutting the axis of abscissas at —B/A. From B/A the value of 


a3/a3 can be found. 


The calculations on which o,4, og are based assume a rigid molecule. The 
effect of the ground state vibration of the molecule was calculated and the 
experimental intensities were modified before comparison with ¢4, og. 


The next section contains the details of the calculation of ¢, and og and the 
subsequent section is the discussion of the vibration correction. 


(c) Cross Sections Averaged over the Thermal Motion 

In the center of mass system, the differential cross section for scattering of 
a neutron of momentum >» through an angle 6 by a deuterium molecule which 
undergoes a transition as the result of the scattering, the rotational quantum 
number changing from J to J’, is (16) 
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oy’ (0, po) = w (2J' + ne 7 (2L + 1) pik (kre), 


where p is the momentum of the neutron after the collision and differs from po 
because of the difference of energies of the rotational states. Assuming a rigid 
molecule having the nuclei at a separation r, 


pm pi+ 4/5 Gar ue =r | + 1) — JI +1) |. 


The values of L are |J — a. : ts — J'| +2,...,J+ J’. The coefficients 
Cz,,,' are discussed in references (8) and (16). k is the magnitude of the vector 
k given by 
k = x(Po — P) 

jx(x) is the spherical Bessel function of order L + 1/2, u is determined by the 
scattering lengths a3/2 and a/2. They enter in the combinations 

(2a3/2+ @1/2)?= A 

(a3/2— @1/2)* = B 
and yp depends on the initial and final states according to Table I. 


TABLE I 
THE QUANTITY #t FOR THE POSSIBLE TYPES OF TRANSITION 


| 
Transition FF x Be 
Ortho-ortho Even - even go (A + (5/4) B) 
Para-para Odd —- odd af (A + (1/2) B) 
Ortho-para Even — odd 5f (3/4) B 
Para-ortho Odd - even ae (3/2) B 


In terms of a function 
Pyy (x) = (2J' + 1) * (2L + 1)Cryy jr?(x) 


the formula may be written 


ory (0, Po) =u . Pyy (kr). (2) 
0 


Formula (2) is fundamental for the subsequent work. It gives the cross 
section in the center of mass system for a neutron of arbitrary momentum pp. 
Transforming to a system in which the molecule is moving with a velocity V 
and the neutron with momentum ,’ equal to the incident neutron momentum 
in the laboratory, and averaging over all V according to the Maxwell distri- 
bution, gives the required cross section corrected for thermal motion. 


Spiers shows how to integrate over two components of the velocity, and the 
resulting formula is 


tie eons [x@ Py» (nen da, (3) 
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where ¢ is the angle of scattering in the laboratory system. The variables are 
expressed in terms‘of the angle a between the vectors po’ and R. If 








E = energy of incident neutrons, 
E, = energy of first rotational state of the molecule, 
T = temperature of the gas, 
then g = * . : 
8 akyl 
y = 25E/4k,T, 


X(a) = sin’a cosec*(a + ¢), 
& = 0.8 cosa — 0.5 sin ¢ cosec (a + 9), 

+ 0.2(E,/E) [(J(J + 1)— J'(J'+ 1)] sin (a + ¢) cosec (¢), 
n = V(E/2E)) sin ¢ cosec (a + ¢). 


To keep all errors in the computations less than the experimental error, it 
was necessary to carry out numerical integrations using Equation (3). For 
each integration 11 values of a were taken at 5 degree intervals, with ao as the 
mid-point where ap is the value of a for which £ = 0. Simpson's rule was applied 
to get the integrals. Although the number of points used in the integrations is 
small, the shape of the curve is such that Simpson’s rule should give the areas 
accurately. Thirty-six partial cross sections were corrected in this way for 
thermal motion. Values for molecules at rest taken from Spiers’ report were 
used for other cross sections which contribute only slightly to the final results. 


Apart from general physical constants the quantities required for the com- 
putation are T the temperature of the gas, E the energy of the incident neu- 
trons, and £, the energy of the first rotational state of the molecule. 


T was assumed to be the temperature of liquid nitrogen, 77.4°K., for reasons 
discussed in Section IIIb below. In any case the temperature has only a small 
influence on the result and a considerable departure from this temperature 
would cause no serious error. 


The value of E used in the calculations was 0.07 ev. Although this differs 
somewhat from the experimental value 0.0724 ev. as measured by crystal 
diffraction, (Section IV), it allowed a direct comparison to be made with Spiers’ 
results which were based on 0.07 ev. Also it was not convenient, at the time the 
computing was begun, to settle on a final value. 


E, was taken to be 0.007416 ev. and determines the internuclear separation 
r, in a rigid molecule by the formula r, = (h/2m)/W MEy = 0.747 X 1078 cm. 
The correction for vibration was based on a molecule having the separation 
ro, the value of r for the maximum density of the zero point vibration. Conse- 
quently ro should have been used in the calculations instead of the r, which is 
implicit in Z;. The relationship is ro>= 1.0056 r-. 





To allow for the above alterations in neutron energy and internuclear separ- 
ation, the theoretical cross sections were taken to apply to angles differing 
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slightly from the angles used in the computing. The basis for this change is that 
kr, is the variable that determines the angular variation of cross section, and 
to a sufficient approximation kr, is proportional to 7 sin (0/2)/A. With the new 
values of r and \, @ was changed to keep this quantity unaltered. The values 
of @ were then transformed to the laboratory system giving the new angles. 


For each transition J — J’ and angle ¢ the quantity computed is o,(¢)/u= 
S;y7. The sum of o,, over J’ is 
(4 + -— a)r Say ot “BY ner even J 
16 4 J’ even 
oJ = — 


+ 5B vo Sin + 58% Syy odd J 
J’ odd J’ even 
as can be seen from Table I. 


The a,’s are added according to the statistical weight of the Jth state. The 
ratio of ortho to paradeuterium was not determined experimentally and con- 
sequently the statistical weights are not known exactly. The extreme cases, 
however, correspond on the one hand to room temperature (2/3 ortho, 1/3 
para) and on the other hand to 77.4°K. (70.3% ortho, 29.7% para). The first 
case will be called ‘‘normal’’ the second case “equilibrium” mixture. The 
fraction of molecules in a state J is either 





fx = Ws for equilibrium mixture 
LW, 
or fy (even) = 2W,/(3 L Ws) ; fz waa) = Wy/(3 & W,) for normal mixture. 
J even J odd 
In all cases W, is computed for a temperature of 77.4°K. 
The cross section per molecule is 
© x oy fy 


and this becomes o@ = Ao, + Bog 


{ ¥ fx Say + a fa x woe” \ 


ae J' even 


{x (Bx Sy + — = ) 
J even 4 De 


+E p35 Sut rin git. 


J’ odd \2 J’ even 


16 
where o4 = 35 
_ 16 
25 


cy:™ 


In Table II are tabulated the S,, for molecules at rest and for a gas temper- 
ature 77.4°K. The percentage contribution of each transition to the 4 and og 
is included. It is interesting to note that at small angles the scattering is almost 
entirely elastic, but that at larger angles inelastic scattering predominates. 
This has the consequence that changes in parity can be of importance only at 
the larger angles. The two sets of o4, og for “‘normal”’ and “equilibrium” mix- 
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tures corrected for thermal motion are in Table III. From the experimental 
results the ratio ef the scattering lengths was found to differ by 5.6% de- 
pending on whether ‘‘normal” or “equilibrium” mixture is assumed. Values 


TABLE II 
THE SCATTERING FACTORS SjJ’ 











| Sys’ Percentage contribution 

+ . Molecules Gas atc 
oO ’ 
at rest 77.4°K. tooa® to oB* 
(a) Scattering angle @ = 10° 
0 | 0 (1. 480)** 1. 496 59.4 68.6 
0 1 (0. 0727) 0.0751 2.1 
0 2 (0.0022) 0.0027 0.1 0.1 
0 3 0. 0005 0.01 
1 0 (0. 0279) 0.0287 0.8, 
1 1 (1. 480) 1.497 29.3 13.5 
1 2 (0. 0583) 0.0613 1.6 
1 3 0. 0062 0.1 0.06 
2 0 0.0005 
2 1 (0. 0429) 0.0445 0.2 
2 2 (1. 480) 1.495 10.6 12.2 
2 3 0.0749 0.4 
2 4 0.0253 0.2 0.2 
3 0 
3 1 0.0018 
3 2 0. 0603 0.02 
3 3 1. 486 0.3 0.1 
3 4 0.0927 0.03 
(b) Scattering angled = 30° 

0 0 (0. 9604) 0.9618 57.2 48.7 
0 1 (0. 4376) 0. 4473 13.6 
0 2 (0. 0344) 0. 0363 2:2 1.8 
0 3 0.0029 0.08 
1 0 (0. 1690) 0.1724 5.2 
1 1 (0. 9770) 0. 9800 28.8 9.8 
1 2 (0. 2762) 0. 2828 8.5 
1 3 0.0240 0.7 0.3 
2 0 0.0136 0.1 0.1 
2 1 (0. 2206) 0. 2254 1.2 
2 ys (0. 9723) 0.9741 10.3 8.7 
2 3 0. 2479 1.3 
2 4 0.0392 0.4 0.4 
3 0 0.0005 
3 1 0. 0226 
3 2 0. 2537 0.07 
3 3 0. 9634 0.3 0.08 
3 4 0. 2256 0.06 


* The Ga, Ox are for equilibrium at 77.4°. 
** Values in brackets ( ) were not used in the calculations. 
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TABLE 11—Continued 
THE SCATTERING FACTORS SyJ’ 








SJ’ Percentage contribution 

- ” Molecul G 

) eh | SL WAR: hs 

(c) Scattering angle = 50° 

| 0 0 (0. 4109) 0. 4054 44.6 24.7 
0 1 (0. 6718) 0. 6803 24.8 
0 2 (0. 1285) 0. 1331 14.6 8.1 
0 3 0.0110 0.4 
1 0 (0. 2542) 0. 2566 9.2 
1 1 (0. 4835) 0. 4807 26.1 5.8 
1 2 (0. 4220) 0. 4280 15.4 
1 3 (0. 0637) 0. 0621 3.4 0.7 
2 0 (0.0511) 0. 0530 1.0 0.6 
2 1 (0. 3286) 0. 3323 2.1 
2 2 (0. 4631) 0. 4595 9.0 5.0 
2 3 0. 3451 2.2 
2 4 0.0519 1.0 0.6 
3 0 0. 0042 
3 1 0.0819 0.04 0.01 
3 2 0. 3671 0.12 
3 3 0. 4596 0.2 0.05 
3 4 0. 3038 0.10 

(d) Scattering angled = 70° 

0 0 (0. 1160) 0.1112 19.5 8.7 
0 1 (0. 6062) 0. 5956 28.1 
0 2 (0. 2212) 0. 2301 40.3 18.1 
0 3 0.0173 8 
1 0 (0. 2211) 0. 2208 10.3 
1 1 (0. 2478) 0. 2470 21.3 3.8 
1 2 (0. 3963) 0.3949 18.4 
1 3 (0. 0933) 0.0931 8.0 1.4 
2 0 (0. 0893) 0. 0920 2.9 1.3 
2 1 (0. 2952) 0. 2953 2.5 
2 2 (0. 2177) 0. 2101 6.5 2.9 
2 3 0. 3204 2.7 
2 4 0. 0408 1.2 0.6 
3 6 0.01421 0. 005 

| 3 1 0. 1364 0.1 0.02 
3 2 0. 3203 0.13 
3 3 0. 2054 0.2 0.03 
3 4 0.11 


0. 2754 | 





* The 0a, Op are for equilibrium at 77.4°. 
** Values in brackets ( ) were not used in the calculations. 


of o,4 and og midway between “normal”’ and ‘‘equilibrium’”’ were finally used, 
introducing at most an error + 2.8%. It is possible that the gas was in such an 
intermediate stage between ‘‘normal”’ and “equilibrium” because at our pres- 
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sure and temperature and in metal containers the para-ortho hydrogen con- 
version takes place in a few days (4, page 61). 


Before computation of the above results was begun an attempt was made to 
avoid the numerical integration of Equation (3). y is so large at 77.4°K. that 
the main contribution to the integral is from the region close to § = 0. If the 
variable is changed from a to é the integrand near § = 0 may be written as a 
power series in £ multiplied by exp (— é). An integration of the first few 
terms between the limits + for £ gives an approximate cross section which 
seemed at first to be a promising alternative to numerical integration. How- 
ever, no simple reliable estimate could be made of the error and this method 
was abandoned, except that cross sections so computed were compared with 
the results of the numerical integrations as a check for gross errors in the latter. 


TABLE III 
THE QUANTITIES 04, OB 























F Equilibrium mixture Normal mixture Final choice 
Scattering angle |__ See eke die TMi ie, coh ce ee 
oA OB OA OB oA ‘OB 
10° 0. 9616 1.041 0.9616 1.017 0. 9616 1.029 
30° 0. 6419 0. 9421 0. 6419 0. 9326 0.6419 0. 9374 
50° 0. 3472 0. 7843 0. 3472 0. 7867 0. 3472 0.7855 
70° 0. 2183 0. 6069 0. 2183 0.6117 0. 2183 0. 6093 


(d) Correction for Molecular Vibration 

The energy difference between vibrational levels is so large (= 0.5 ev.) that 
only the ground state occurs in the gas, and the neutrons have insufficient 
energy to cause transitions to higher vibrational levels. The scattering may 
then be treated as being off the probability density of the ground state vibra- 
tion. This interpretation is in accord with the result of retaining the vibrational 
wave functions in Equations (10) to (20) of reference (8). (An analogous case 
is the elastic scattering from a single rotational level as may be seen from 
Equation (19) of reference (8).) 

For a Morse potential (13, page 273) the radial probability density of the 
ground state vibration is given by 


dp me eek ll ae :) 
Y dr 4/2 (1+ #") exp ( = dr, (4) 
where u = a(r—?o), 
ro = separation of maximum density, 
b = 46.56, 
e= 1.49/ro. 


The constants of the deuterium molecule were obtained from reference (4) 
Table 54. 

The difference between the scattering from diatomic molecules having line 
scatterers with the radial distribution of Equation (4), and diatomic molecules 
having point scatterers at a fixed separation 79, was calculated with the result 
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0.00244 10.333 x? + (0.98 + x?) Si) * + 0.98 cos x} 
x 


Sesh 


* tare 


Us = 


(5) 


x 
Ix is the intensity scattered by the vibrating molecule, 
Io is the intensity scattered by the rigid molecule, 
eo 4x ro sin 0/2 
5\/4 
6 = angle of scattering in the center of mass system, 
5\/4 = neutron wave length in the center of mass system. 


The upper signs refer to scattering without change of parity, the lower signs 
to scattering with change of parity. A semiclassical approximation was used 
instead of an accurate summation over rotational transitions. This is the same 
as the method used for calculating X-ray scattering (2, page 159). The mole- 
cule was treated as a rigid structure of scatterers, and the intensity of the wave 
scattered in the direction @ was calculated and then averaged over all 
orientations of the molecule. The calculation was done in the center of mass 
system and subsequently transformed to the laboratory system. This method 
was justified experimentally as an approximation for molecules without nuclear 
spin (1). Spiers (16, 17) has given a theoretical discussion and has shown that 
with deuterium, a distinction is necessary between transitions in which the 
parity is unchanged and transitions in which the parity changes. In the former 
case the two nuclei scatter identically; in the latter case, a phase difference of 
180° is introduced. 


The fraction U by which the scattered intensity differs from that calculated 
for a rigid rotator is the weighted mean of U, and U_. The semiclassical formu- 


lation (16) for a rigid rotator gives (1/2) (1 + se) = sum of 7, over J’ of 
x 


same parity as J, (1/2)y (1 - =) = sum of o,, over J’ of different parity 
x 


from J. With the obvious notation yoo= yu for ortho-ortho transitions, etc., and 


with wt = (1/2) (1 + sx) , the differential cross section is 


Le fs & esr = aw, + Bw, 
where a = woo > fy + upp Dd fy, 


o 


J even J odd 
B = pop D fa + upo X fi. 
J even J odd 


The effect of vibration is to change o to o(1 + U) and ws to wi(1 + Ux). As 
intensity is proportional to o it follows that 


ki aw,U, + Bw U- 
aw, + Bw- ; 
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The values of Uz and U are plotted in Fig. 1 as functions of the angle of 
scattering in the laboratory system. 





SCATTERING ANGLE 


Fic. 1. The effect of ground state vibration of the molecules. The three curves show the change 
required in the scattering calculated for a rigid rotator model to yield the correct scattering. Ut applies 


when there is no change in parity, U_ when there is change in parity, and U is the average for the 
gas at 77.4°K. 


It is possible to avoid the semiclassical calculation and correct the individual 
o,y7’s from a comparison of j,(kr.) with the average of jz (kr) over the vibration. 
The increased labor was ‘not justifiable except for small angles (< 30°) where 
a good approximation is readily obtained. In this region close agreement was 
found with the semiclassical values. This gives confidence in ia 


For convenience, the experimental results were modified for comparison with 
the rigid rotator model, i.e., the experimental results were divided by 1 + U. 


III. Apparatus 
(a) Spectrometer 


Fig. 2 shows the experimental arrangements used for the gas scattering 
experiments. A beam of neutrons from the reactor impinged on a sodium 
chloride crystal set to diffract neutrons having a wave length of 1.063A. The 
diffracted beam passed through the scattering cell where a fraction was scat- 
tered by the gas. A boron trifluoride counter was used to measure the scattered 
intensity as a function of scattering angle over the range 11° to 67°. Steel 
collimators defined the primary, diffracted, and scattered beams. They were 
81 in., 13 in., and 14 in. long, respectively, and had 1/2 in. by 1 in. rectangular 
holes cut down their lengths. Surrounding the first and part of the second 
collimator was a fast neutron and y-ray shield of paraffin and barytes and a 
further layer of steel. Surrounding the third collimator and the counter was a 
neutron shield of 8 in. of paraffin and 1 in. of boron carbide. 





HURST AND ALCOCK: THE SCATTERING LENGTHS OF THE DEUTERON 47 


The first collimator formed part of the “hole extender”, a 12 ft. tube used 
to offset the spectrometer from the congested area of the reactor face. Internal 
details of this extender, not shown in Fig. 2, were irises of boron carbide and 
lead to reduce radiation in the neighborhood of the extender, and a 12 ft. 
evacuated tube to prevent reduction of the beam intensity by air scattering. 













EXPERIMENTAL HOLE EXTENDER 
STEEL COLLIMATOR 


PARAFFIN, BORIC ACID NEUTRON BEAM FROM PILE 


@ BARYTES MIXTURE /| Px 7 


NaCl CRYSTAL (10,0) PLANES 


NEUTRON BEAM OF WAVELENGTH |+-063A 







STEEL COLLIMATORS BF, COUNTER 


SCATTERED NEUTRON BEAM 


PARAFFIN 
Fic. 2. Neutron crystal spectrometer as arranged for gas scattering. 


The neutron counter, an enriched boron trifluoride proportional counter 
2.5 in. diameter by 25 in. long, has been described previously (7). It should be 
noted that the counting efficiency is almost independent of neutron energy in 
the region 0.07 ev. This characteristic was important because the energy of a 
scattered neutron depended on the scattering angle, the translational energy 
of the scattering molecule, and any change in rotational energy of the molecule. 
The range of neutron energies, however, was such that differences in counting 
efficiency could be ignored. 


The spectrometer was operated automatically. A servo-mechanism drove the 
counter arm through the chosen angular range, setting it at regular 8° intervals 
with a repeatable angular accuracy of + 2’. The servo-mechanism was trig- 
gered by the last scaler of a neutron counting circuit, after a preset number of 
counts. The time taken at each point depended on the reactor power level 
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t 
during the counting interval. Therefore | Ndt was recorded, where N was 
0 


proportional to the neutron flux in the primary beam. This was done by re- 
cording for each point the total number of counts in a neutron counter located 
vertically above the spectrometer crystal. This monitor provided adequately 
for small fluctuations in reactor power. To ensure, moreover, that the measure- 
ments were made under the steadiest conditions, counts were allowed to accu- 
mulate only when the reactor was on automatic control. 


Cu TUBING TO PRESSURE GAUGE, 
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Fic. 3. Scattering chamber. 


(b) Scattering Chamber 

‘The scattering chamber shown in Fig. 3 was a double-walled steel Dewar 
flask with an upper compartment for the liquid nitrogen, and a lower cell for 
the deuterium gas. This gas cell was 2 in. in diameter and 2.5 in. high with a 
0.030 + 0.0005 in. thick polished steel wall. Internal cadmium-coated steel 
baffles served two purposes. Primarily they reduced double scattering, but also, 
since they were firmly screwed to the underside of the liquid nitrogen com- 
partment, the baffles helped maintain the gas at a uniform temperature of 
77.4°K. The Dewar was continuously evacuated with an oil diffusion pump to 
about 3.10~4 mm. of mercury pressure. 


At small angles (< 30°) the neutron counter was exposed to part of the wall 
of the gas cell which was in the neutron beam, (Fig. 3, Section A-A). This, of 
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course, increased the background of slow neutrons. However, the effect was not 
as great as might be expected, for the steel wall was crystalline and most of its 
scattering, like that of a crystalline powder, was into Debye—Scherrer cones. 
No cone was inside an angle of 30°, and there was little diffuse scattering. 


An Ashcroft Laboratory test gauge, 0-1000 p.s.i. range with 5 p.s.i. sub- 
divisions, and evacuating and filling needle valves, were connected to the gas 
cell with a length of copper tubing. The pressure gauge was calibrated by the 
Design and Development Branch at Chalk River and found to be accurate. 
Liquid nitrogen was blown into the Dewar from a 25 liter flask at 15 min. inter- 
vals by means of an electrically actuated compressed air valve, the valve being 
opened by a time switch and shut off by a tungsten coil (Fig. 3) connected into 
a balanced bridge. As a precaution against failure of the liquid nitrogen system, 
which would produce a dangerously high gas pressure, there was in the system 
a rupture disk rated at 525 p.s.i. maximum. An external heater prevented dew 
formation on the outside of the Dewar at the gas level. 


IV. Experimental Results 


(a) Angular Distributions 


_ The intensity of neutrons scattered from deuterium gas was measured at 
eight equally spaced angles: 11.2°, 19.2°. . . . 67.2°. One thousand counts were 
recorded at each point and the whole angular range was swept back and forth 
through 16 half cycles, the whole run lasting 64 hr. 


Background runs, with the gas cell evacuated and at room temperature, 
were made before and after the gas scattering run. The first background run 
was five half cycles, the second three, with 100 counts at each point. 


A separate measurement was made of the fast neutron background alone, by 
substituting a composite absorber of sheet cadmium, paraffin, and boron car- 
bide for the gas cell. With 100 counts at each point, a run of three half cycles 
was made. 


To calibrate the angular scale, a transit was mounted on the axis of rotation 
of the counter arm. The telescope was aligned on a distant point, and the 
counter arm and transit scales were correlated through an angle of 90°. Zero 
angle was determined by measuring the shape of the monokinetic beam, and 
finding the mid-point. A neutron absorber was used to reduce the intensity of 
the main beam. This was done both before and after the angular distribution 
measurements, and the two measurements agreed. 


The measured angular distribution had a 12:1 variation in counting rate. 
This large variation was due to the true deuterium pattern, the geometry of 
the gas cell and collimator slit system, and the background. The total back- 
ground varied 30: 1; but the fast neutron component was constant. The count- 
ing rates are illustrated in Table IV for three angles: 11.2°, 35.2°, 67.2°. These 
rates are expressed as counts per minute for clarity, though in all the calcula- 
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tions counts per 10‘ monitor counts were used. The gas was at an average 21 
atm. pressure and 77.4°K. The errors quoted are standard deviations. 


TABLE IV 


MEASURED NEUTRON INTENSITIES SCATTERED FROM DEUTERIUM GAS 
(21 ATM. AND 77.4°K.) AS A FUNCTION OF SCATTERING ANGLE 


Scattering angle 


Gas-scattered + Total background (c/m) 
Total background (c/m) 
Fast neutron background (c/m) 





(6) Wave Length Measurement 


The neutron energy in the monokinetic beam from the crystal spectrometer 
was measured before and after the gas scattering experiments. For this the gas 
cell was replaced with a sodium chloride crystal, and rocking curves were made 
on both sides of the beam from the first crystal. The narrow rocking curve 
obtained when the planes of the two crystals are parallel is a measure of the 
imperfections of the crystals, while the broad curve obtained when the crystals 
are antiparallel is a measure of the energy spectrum of the diffracted neutron 
beam. The angle between the two peaks is four times the Bragg angle 03, and 
the neutron wave length is \ = d sin 0s, where d = 5.628A is the lattice 
spacing of the crystal. 


The angle between the maxima of the diffracted peaks was 43.56°, corres- 
ponding to a neutron wave length of 1.063A or an energy of 0.0724 ev. The 
width at half maximum of the broad diffraction peak was 1.15° and of the 
narrow peak 0.22°. Therefore, the spread in wave length of the diffracted 
neutron beam from the spectrometer was approximately + 0.05A. 


(c) Gas Purity 

Deuterium gas, drawn from the gas cell at the end of the angular distribution 
run, was analyzed for isotopic composition by the Chemical Control Labor- 
atory at Chalk River, using a Nier-type 60° mass spectrometer. The average 
result of four samples taken was Dz molecules 99.20%, HD molecules 0.80 + 
0.06%. 


(d) Transmission Measurements 

Transmission measurements were made on the scattering gas and cell walls 
as a rough over-all check on these materials. For these measurements, the 
counter arm was set to 0°, and the beam was stopped down with a cadmium 
iris. 


The scattering cross section per atom of deuterium gas at 21 atm. pressure 
and 77.4°K. for 0.0724 ev. neutrons was found to be 3.4 + 0.4 barns. This is 
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in agreement with the accurate cross section for neutrons of a few volts energy: 
3.44 + 0.06 barns (6). The scattering cross section of iron for 0.0724 ev. 
neutrons was found to be 12.6 + 0.5 barns, in agreement with the known value 
12.7 barns. 

V. Reduction of Experimental Results 


Before the measured angular distribution can be compared with theory, 
attention must be given to background, gas pressure, effective scattering 
volume, double scattering, and hydrogen impurity. These complications are 
treated below in the order of their application. 


(a) Background 

Background arose from two sources: fast neutrons which penetrated the 
counter shield; and slow neutrons which had been scattered into the counter 
collimator by the steel walls of the gas cell. The fast neutron component was 
unaffected by the addition of gas to the cell. The slow neutron component was 
assumed to traverse the cell, and the measured value was multiplied by exp 
(— ds), where d is the cell diameter and s is the fraction of beam scattered per 
unit length by the gas. This correction for gas transmission is appreciable at 
small scattering angles but is only 0.5% of the gas-scattered intensity for angles 
greater than 30°. After this adjustment the background was subtracted from 


' the total intensity to give the gas-scattered intensity. 


(b) Gas Pressure 

A small leak in the system caused the gas pressure to drop steadily through- 
out the experiment, from 22.5 to 19.5 atm. Frequent readings of gas pressure 
were made and a smooth curve of pressure versus time has been drawn from 
which the exact gas pressure for each 1000 count interval was determined. The 
gas-scattered intensity for each interval was then corrected for any pressure 
deviation from an average 21 atm. 

The 15 min. liquid nitrogen filling cvcle -aused a varying temperature 
gradient along the filling pipe and therefore a fluctuation in the gas pressure 
of about plus or minus 1‘%. Since the average error resulting from this ripple 
must be very much less than 1%, it has been neglected. 


(c) Scattering Volume 
The effective scattering volume of gas as a function of the scattering angle is 


a valid) 


2r(v - 2) + —————_ | for 0 < @ X 2 
2 $ 


cos — 
| 
w T 


os for 2y < @< — 
sin 2 





where y = sin}. V is the volume of gas of unit height, 7 is the radius of 


7 
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the gas cell, w is the collimator slit width, and ¢ is the scattering angle. The 
gas-scattered intensities were therefore divided by V to give the scattering per 


unit volume. 


(d) Double Scattering 

An approximate correction for double scattering in the gas, and between the 
gas and cell wall, was made. The correction was based on two very crude 
assumptions: that the gas scatters isotropically in the laboratory system, and 
the wall scatters not at all for scattering angles < 30°, isotropically for 
angles > 30°. The necessary integration was made numerically for the 
geometry of the experiment, i.e., 1/2 in. collimated beams, 2 in. I.D. gas cell, 
and internal nonreflecting baffles as shown in Fig. 3 (Section A-A). The result 
may be written: 

Ip(¢) = [aP(¢) + 5Q(¢)] Zs, 

where Jp is the intensity due to double scattering as a function of the scattering 
angle ¢; J, is the average intensity due to single scattering; a is the fraction 
of beam scattered by the gas, at the temperature and pressure of the experi- 
ment; } is the fraction of beam scattered by one wall; P and Q are functions 
of ¢ as given in Table V for the scattering angles 0°, 30°, 45°, 60°, and 90°. 





TABLE V 
DouBLE SCATTERING TERMS 
Scattering angle (@) P(o) 0?) 
0° 0.177 0. 229 
30° 0. 235 0.114 
45° 0. 269 0. 108 
60° 0.305 0.099 
90° 0. 280 0.073 


Ip(¢)/I, was found to be nearly independent of angle for the conditions of 
this experiment, varying between 3.1% and 3.4% over the angular range 11° 
to 67°. I, corresponds approximately to the experimental intensity at 67.2°. 
This was determined by a numerical integration over 4x steradians of Spiers’ 
differential cross section ¢4 + ; OB; 7 being found by a successive approxi- 
mation. Finally Ip was subtracted from the scattered intensity per unit volume, 
leaving an experimental intensity I(¢). 


(e) Hydrogen Impurity 

The scattering cross section of hydrogen is approximately six times that of 
deuterium, and since the deuterium gas of this experiment contained 0.4% 
hydrogen, a correction for this impurity was required. 

At such a low concentration the hydrogen would be entirely in the form of 
HD molecules. Dr. N. K. Pope of this laboratory has kindly calculated the 
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differential cross section for the HD molecule and from his general formula has 
derived the semiclassical approximation, viz. 


a, (8) = “4 (4% + Ap + 0.75 By + 2Bp + 2AyAp sin) (6) 





4nr,.sin . 
2, 
where x= — 
4/3 
re = 0.747A, 
@ = angle of scattering in center of mass system, 


\ = neutron wave length in laboratory system = 1.063A, 
Ay = 3a;+ ao= — 0.75 X 107" cm., 
By = 2(a;— ao) = 5.81 x 10-2 cm., 
Ap = 2a3/2+ a1/2= 1.346 & 107” cm. 
Bp = 43/2— 41/2= + 0.567 X 10~™" cm. 
a, and ap are the scattering lengths for hydrogen and are given most accurately 
by the total scattering cross section oo, and the coherent scattering amplitude 


ay(10, 12). 
oo = 42(3/4 ay?+ 1/4 ac?) = 20.36 & 107% cm?. 
Qy = 2(3/4a, + 1/4 a0) = — 0.375 X ~"? cm. 


3/2 and ai/2 are the scattering lengths for deuterium given in this report 
(Section VI). They were applied in a successive approximation. 


The measured intensities are the sum of the intensities scattered by the HD 
and Dz molecules. If J(@) is the experimental intensity as given in (d) above, 
and I(Dz, ¢) the intensity due to deuterium molecules. 
1 
I(D2, ¢) = I 
(Dz, ¢) (¢) i+ Ng 
where g = oxp(¢)/cp2 (¢), 





y = 0.008 

0.992 
op2(¢) are the cross sections for scattering in the laboratory through an angle 
¢ by HD and D; respectively. The semiclassical formula was used not only for 
op, but also for ope, since the ratio of the two would then minimize differences 
between the approximate and exact formulas. The HD cross section of Equa- 
tion (6) is independent of J and we have only to transform ¢,(@) from the 
center of mass to the laboratory to get ¢yp(¢). The deuterium cross section 
had to be summed as described for the vibration (Section II(d)) and then 
transformed to the laboratory system. The correction increases smoothly from 
2.6% at 11° to 5.2% at 67°. 


, is the number of HD molecules per D2 molecule, andeyp (¢), 


(f) Outer Effect 
Owing to the high pressure of the gas, the scattering at small angles is less 
than that for free molecules. Debye (3, 14) calls this the “outer” effect and 
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attributes it to the unavailability of the space occupied by the molecules. For 
X-ray scattering by a gas of diatomic molecules which act like spheres of 
diameter 2 in collisions, but like pairs of particles separated by a distance r, 
in scatter..ig, Deb: has shown that the intensity is proportional to 


, - sa se 
(1/2) (4+ sinx\ 4nD*N sin x/2\?sin u -- u cos u (7) 
x x us 
2 
4nr, . 0 
where. < = ——‘gn— , 
r 2 
u = Dx/s,, 
N = number of molecules per unit volume, 
>’ = wave length, 


6 = scattering angle. 


This formula applies also to neutron scattering by deuterium if transitions 
involving parity change are unimportant and the semiclassical approximation 
is assumed. Both these conditions hold at small angles (< 20°). 


The ratio of the second term to the first was computed, using D = 2.72A 
(11, page 183). The computation was made for the center of mass system, and 
then transformed to the laboratory system. The ratio was then applied to the 
calculated intensity with the result shown in Fig. 6, to give much better agree- 
ment with the experiment. It should be noted, however, that neither of the 
points affected by the ‘‘outer”’ effect was used in the plot from which a;/2/a3/2 
was found. 


VI. Determination of Scattering Lengths 


(a) Ratio of a1/2 to a3/2 

The final corrected experimental intensities, modified also for molecular 
vibration (Section II(d)), are given in Table VI. In Fig. 4 these experimental 
intensities, divided by og, are plotted against o4/og (Section II(b)) for the 
angles 27.2, 35.2, 43.2, 51.2, 59.2, and 67.2 degrees. The ordinate scale is arbi- 
trary and is therefore not numbered. Two abscissa scales are shown, o,4/a, and 
the corresponding scattering angles. The most probable straight line has been 
drawn through the points, and cuts the x-axis at —B/A = — 0.173. A stan- 
dard deviation has been assigned and is shown by two further straight lines. 
B/A is then 0.173 with a standard deviation of + 0.022 and a probable error of 


TABLE VI 


EXPERIMENTAL SCATTERED NEUTRON INTENSITIES 
(with all corrections applied) as a function of scattering angle 


—— —= 
27.2° | 35.2° | 43.2° | 51.2° | 59.2° | 67. 2° 





Scattering angle | 11.2" 19. 2° 











586.3) 578.2; 479.2 392.4| 316.0) 255.6) 215.4 185.3 
+ 12.9)+ 9.2 5.7)+ = 3.2)4 2.6) ae 1.9 
' a { | i 


| 


Intensities with standard 
deviations 
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Fic. 4. Plot of experimental intensities — oB versus OA/OB. 























SCATTERING ANGLE in LABORATORY 


Fic. 5. Theoretical differential cross sections corrected for thermal motion of molecules. 
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+ 0.016. To the statistical error, + 0.016, must be added errors due to uncer- 
tainties in the theory or experimental conditions; in particular the double 
scattering correction is inexact and the ortho-para ratio is unknown. Finally: 


2 2 
B/A = 0.173 aoe Since B/A = (1 ~ 01/»/012) /(2 + a1/»/03/) , 


it follows that: @1/2/a3/2= 3.2 + 0.3 or 0.12 + 0.04. 












23 Twene : ors. “— 





SCATTERED NEUTRON INTENSITY 
§ & 
ob 


SCATTERING ANGLE in LABORATORY SYSTEM 


Fic. 6. Angular distribution of scattered neutrons from deuterium gas at 77.4°K. and 21 atm. 
pressure. 


Fig. 5 shows the calculated differential cross sections o4 and og, together 
with the cross section ¢,-+ 0.173 og, where all three curves have been normal- 
ized to unity at 0°. The curves are based on the accurate values of Table II. 
Spiers’ results for molecules at rest (but recalculated for 77.4°K) were used 
for interpolation. Fig. 6 shows the same ¢4+ 0.173 og, as a full line curve, with 
the experimental intensities of Table VI fitted at 27.2°. The full line curve is 
applicable for low pressures only. According to formula (7), it is changed at 
low angles to the dotted line for the average operating pressure of 21 atm. 
but is otherwise unaffected. 


The ratios of the scattering lengths were combined with the cross section of 
the free deuteron as measured by Fermi and Marshall (6), viz., 3.44 + 0.06 
barns, to give the scattering lengths. The knowledge that the coherent scatter- 
ing length is positive (15) reduced the four possible values of each length to 
two values which are 
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@1/2 = 0.826 + 0.012) 0.07 + 0.03 | 
oe hee ~ < 107% cm. 


a3/2 = 0.26 + 0.02 | 0.638 + 0.006 


(b) Previous Measurements 


The methods for the determination of scattering lengths are compared on 
Fig. 7. With a1/2 and @3/2 as co-ordinates a value of o7 establishes an ellipse. 
A ratio @1/2/a3/2 is a straight line and the intersections with the ellipse give the 
scattering lengths. Experimental uncertainties spread the lines into bands and 


a; [y= 0.2 +0.04 


Oo (D,) = 21.3 to2 





-0. ° Ol 0.2 0.3 04 05 06 0.7 os os LO 
OQ; 116*om,) 


Fic. 7. Comparison of methods for determination of the scattering lengths. The total cross 
section or gives an ellipse and the scattering lengths are obtained from the intersections of this 
ellipse with another experimental quantity. Three such quantities are (i) o(D2) the cross section 
of deuterium gas for low energy neutrons, Fermi and Marshall (6), (it) O¢ the coherent scattering 
cross section, Wollan, et al. (18), (tit) a, a3 from the present experiment. 


the intersections become areas. The analysis of the experiment of Fermi and 
Marshall (6) is equivalent to finding the intersection of 


o(D2) = 6.271 (2a3/2?+ @1/2")— 1.865 (a3/2 — a1/2)* \ 


with or. This is shown on the same figure. Crystal diffraction (15, 18) gives the 
coherent scattering length fp (usually quoted in the form of a cross section 
o.= 4xfp) from which the straight line 2a3/2+ @1/2= 2fp is obtained. 


In principle the intersection of any two experimental lines on Fig. 7 yields 
a possible pair of scattering lengths, but the total cross section, being the 
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simplest to measure, is taken as fundamental. The shaded areas mark the 
regions of experimental values of a3/2 and a1/2 according to the three investi- 
gations, all referred to o7= 3.44 + 0.06. It will be seen that the values of a3/2 
and 4;/2 from gas scattering lie well outside the limits set by Fermi and Mar- 


shall. They lead to a coherent scattering cross section of 5.8 + 0.2 barns, and 


so are consistent with those deduced from the Oak Ridge value of 5.2 + = 


barns; but they are more closely defined. 
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THE PRECISION MEASUREMENT OF HALF LIVES OF 
RADIOACTIVE SUBSTANCES! 


By Lorna M. SILVER? 


Abstract 


A vibrating reed electrometer with feedback amplifier has been used in con- 
junction with a differential ionization chamber to measure accurately the decay 
periods of several radioactive isotopes. This method proved to be simple, 
reliable, and free from uncertain corrections, with a low background due to the 
differential technique adopted and a low drift rate (normally less than 5 X 107!7 
amp.) due to the feedback electrometer circuit. This high and easily varied 
sensitivity enabled the decay to be followed over 12 periods, making possible the 
detection of, and correction for, any weak long lived impurity activity, thus 
eliminating one major source of error in half life measurements. Least squares 
analysis of the data showed that a quadratic mean error as small as 1/2% can 
be easily achieved with this apparatus. Half life values were found as follows: 


a9Cu® 12.88 + 0.03 hr. 
2Cu®é 4.34 + 0.03 min. 
2g Ni® 2.564 + 0.005 hr. 
7Au' . 2.69 + 0.01 days 
agIn'6 53.93 + 0.03 min. 
Introduction 


It is surprising to find that only 10 isotopes in the most recent Seaborg— 
Perlman Tables (12) have their half lives quoted to four significant figures, 
even though an accurate knowledge of decay periods is useful, for example, for 
radioactive standards, or in the monitoring of neutron fluxes and, more gener- 
ally, for sorting out the various isotopes produced in an irradiation. This lack 
of precision in the present values of half lives is due partly to the methods of 
measurement usually employed. Techniques using Geiger—-Miiller counters are 
inadequate since the maximum and minimum rates of counting that are pos- 
sible enable the decay to be followed for only six or seven periods and hence 
make it difficult to detect the presence of a weak long lived activity. Also the 
necessary corrections for dead time are a little uncertain, and the behavior of 
the counters is apt to vary during a run of several weeks. Peierls (10) has 
described the best method of taking counter data and subsequently evaluating 
the mean life of a substance, but it requires continuous recording and constant 
counter behavior during the run, so has not often been adopted. Ionization 
chambers with mechanical quartz fiber electrometers are more frequently used 
as they are more stable and no dead time correction arises. In their most usual 
form these are limited by the range of activities measurable and may suffer 
from mechanical zero drift, while highly sensitive instruments are difficult to 
set up. Vacuum tube electrometers for measuring ionization currents can be 
made quite sensitive and capable of covering a wide range but suffer from zero 
drift and for accuracy require continual calibration checking during a measure- 

1 Manuscript received August 23, 1950. 
Contribution of the Chalk River Laboratories of the National Research Council of Canada. 


The experimental work was completed in August, 1948. Issued as N.R.C. No. 2288. 
2 Now at the University of British Columbia, Vancouver, B.C. 
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ment. Indeed, the early determinations of the period of Rn” by Bothe (2) and 
by Curie (3) still represent two of the most accurate measurements ever made. 
The presence of a weak long lived activity is often the main source of error, 
and to detect and correct for this it is essential to use a technique that com- 
bines an extended range with a high sensitivity. 


Experimental Method 


The apparatus used in this work consisted of a differential ionization chamber 
and a vibrating reed condenser electrometer with feedback amplifier of the type 
described by Palevsky, Swank, and Grenchik (9). This electrometer is a very 
stable, wide range, sensitive instrument, its steady drift rate normally being 
less than 5.10~" amp. A differential ionization chamber has the virtue that by 
careful balancing of the capacities of the two halves and thorough cleaning, the 
differential ionization current due to background activity can be practically 
neutralized. A scale diagram of the chamber is shown in Fig. 1. It is open to 
the air, thus necessitating a check on pressure, temperature, and humidity 
variations. The slots for the source holders enabled specimens to be replaced 
with exactly reproducible geometry, as was substantiated by the linearity of 
the results for a Zn® monitor source. The electrometer head containing the 
vibrating reed and shorting switch was mounted on the heavy rigid steel case 
of the chamber. The experiment was carried out in an underground, air con- 
ditioned room where variations in temperature, humidity, and pressure were 
found to be insignificantly small when checked by following the decay of the 
relatively long lived Zn® source. 


Measurements of the ionization current were made by following the rise in 
potential of the chamber collection electrode with the vibrating reed electro- 
meter by timing the pointer movement of a millivoltmeter having a 6 in. mirror 
scale, with a calibrated stop watch reading to 1/100 sec. Before using, the 
amplifier was carefully checked for linearity of response on the various ranges 
but the confirmatory test was the linearity of the results for a simply decaying 
isotope, In''6, It proved possible to follow the decay over 12 half lives, corres- 
ponding to a maximum range of 8000 mv. per min.—limited by timing ac- 
curacy—down to the background activity level of 0.8 mv. per min. Below 0.8 
mv. per min., statistical errors due to fluctuations in ionization current become 
comparatively too large. The electrometer was operated continuously over the 
complete period of several weeks during which the measurements were taken 
in order to avoid errors due to drifts while warming up. 


Sources were prepared by the irradiation of chemically pure samples in the 
Chalk River pile. In cases such as Cu* and Cu® the neutron flux and time of 
irradiation were adjusted to give the maximum ratio of desired to nuisance 
isotope abundances. The maximum initial activity corresponding to the 8000 
mv. per min. deflection of the millivoltmeter was of the order of magnitude of 


10 microcuries. 
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Fic. 1. Differential ionization chamber. A—Connection to electrometer. B, C—H.T. con 
nectors. D—Specimen holders. E—Ring collecting electrode. F—H.T. cylinder. G—Short- 
ing key. 


Results 


The decay periods were obtained by a least squares analysis of the data. 
Table I demonstrates this method of analysis in finding the most probable 
values of y and 4 satisfying the decay equation: 


log (AN) = y — A4, 


where y = constant X log (ANo) 
’ = 0.43429 vA 
X = decay constant = 0.69315 os 0.69315 - 


half life T 
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The value of the rate of charge at any particular time represents the average 
of several contiguous readings taken around that time. Where necessary this 
average was corrected for the finite time interval At, over which the batch of 
readings was taken by the usual equation: Experimental value of activity (at 
interval’s midtime) divided by corrected value of activity (at interval’s mid- 


time) = 1 + (Aas)?/24. 


The table gives for the most probable values 


y = 3.9140 
= 0.023363 
T = 12.88 hr. 


From the residuals, the most probable value of the quadratic mean error in 
Nr Gr BIVES TF = (0.43429) (0.69315) (1 + @/2’)/2’ 
= 12.88 + 0.03 hr. 


TABLE I 
SPECIMEN ANALYSIS OF DATA—CU“ 


Rate of charge,* Midtime, 

millivolts/min, Log (rate) hr. Equations of condition Residuals 
(Proportional to \N) 

7712. 88 3. 8872 1.3178 | 3.8872 = y — 1.3178’ | —.004009 
2599. 22 3.4148 21.597 3.4148 = y — 21.597’ | —.005425 
762. 41 2. 8822 44.071 ete. +.002164 
235. 40 2.3718 65.911 +.002309 
122. 28 2.0874 77. 757 +.009994 
57. 233 1.7576 92.021 +.006447 
46. 172 1.6644 95. 922 +. 008566 
32. 556 1.5124 102. 936 — .003260 
16. 642 1.2212 115. 661 — .009404 
12. 356 1.0919 121. 109 — .007378 


* = Corrected for background of 0.8 millivolt per min. 


In the analysis of Cu“ (Table I), the change in sign of the residuals shows 
that the relationship between time and log (activity) is not linear, indicating 
the presence of a weak long lived activity—perhaps from some impurity in the 
copper sample. The residuals are much more sensitive to an impurity activity 
than any graphical plot of the relationship would be, as can be seen for example 
by comparing Fig. 2 and Table II, where the results for Ni® are shown. In this 
case presumably Ni® has been formed as well as Ni®. The critical deviation 
of the final point was calculated from the statistical fluctuations of the readings 
at this near-background activity level, and the table shows how sensitive the 
precision is to the amount of long lived isotope assumed present. All points 
have been assigned equal weights in these analyses. 


A list of the values obtained for the half lives of a few isotopes is given in 
Table III together with values obtained by other workers. 
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Fic. 2. Correction for long lived activity in Ni® sample. 
TABLE II 
DETERMINATION OF BEST FITTING CURVE IN FIG. 2 
| 
Points used Activity (rate of charge) Half life % Error 
A, 5 €, BEF Total, i.e., neglecting the correction 2. 730 hr. 2.81 


for 36 hr. activity 





ABCD: Ee Total 2.628 0.73 


A, B, C, D, E, F, G | Allowance made for 36 hr. 
activity assumed to correspond to 
measured value at time G, and 
activity taken as: 











(1) Mean of experimental rates—| 2.614 0.83 
critical deviation 
(2) Mean 2. 564 0. 192 
(3) Mean + critical deviation 2.500 1.38 
Best value 2.564 + 0.0065 hr. 
Conclusions 


The results indicate that it is possible to measure half lives quickly and 
easily to 1/2%. For the greatest precision after the singleness of the decay has 
been ascertained, one could use the double source differential chamber tech- 
nique (3) in which, after establishing the small correction factors needed for 
pressure, temperature, and humidity, the precision is then limited only by the 
accuracy of the timing device and the strength of the available sources. The 
method could also be used to advantage for measurement of long periods and 
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in experiments to detect minute changes in half lives through changes in elec- 
tron density around the decaying nucleus. Here the ionization current from the 
source in one compartment is balanced against a ‘constant’ activity (such as 
radium) in the other and the unbalanced activity measured at a later time, as 
recently described by Segre and Wiegand (13). For ordinary purposes, however, 
the simple experimental procedure described here enables a decay period to be 
measured to 1/2% with some assurance that there is no error due to a weak 
impurity and without any uncertain corrections. 


TABLE III 





PRECISION VALUES OF SOME HALF LIVES 


























Isotope | Half life found Other values Method Author 
in present work 
29Cu* 12.88 + 0.03 hr. | 12.62 + 0.1 hr. Ionization chamber} A. G. Ward* 
12.8 + 0.1 Pressure chamber| Van Voorhis (14) 
with F.P. 54 am- 
plifier 
aCu®® | 4.34 + 0.03 min. | 5 min. Pressure chamber| Amaldi, others (1) 
and Perucca elec- 
trometer 
wNi® | 2.564 + 0.005 hr.) 2.6 + 0.03 hr. Livingood and 
Seaborg (7) 
voAu™ | 2.69 + 0.01 days| 2.72 + 0.05 days| Ionization chamber] A. G. Ward* 
2.73 + 0.02 Geiger—-Miiller Diemer and 
counters Groendijk (4) 
2.7 Electroscope McMillan, Kamen, 
Ruben (8) 
2. 69 D. Saxon (12) 
ain"6 |53.93 + 0.13 min. | 54 min Amaldi, others (1) 
54 Lawson, Cork (6) 
54.31 Russel and Watten- 
berg (11) 
54.05 + 0.16 Geiger-Miiller Graves and 
counters, method} Walker (5) 
independent of 
counting losses 
54.4 + 0.2 Ionization chamber} A. G. Ward* 
54.5 + 0.3 Geiger counters un-| B. W. Sargent** 
corrected for small 
losses 
54.02 + 0.03 Bobin and Lockett 


* Personal communication from Mr. A. G. Ward, Atomic Energy Research Establishment, 


Harwell, England. 


** Personal communication from Dr. B. W. Sargent, Chalk River Labs., NRC, Canada. 
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NOTE ON THE RELATION BETWEEN FEYNMAN’S 
FORMULATION OF SCATTERING PROBLEMS 
AND THE BORN APPROXIMATION! 


By J. R. H. DEMPsTER 


Abstract 
Feynman's formulation of scattering blems is des.) ed briefly, and is 
compared with Born’s approximation. e two proc’/».es are shown to be 
equivalent and to be valid under the same conditic« Feynman’s formulas 


are used to calculate the first and second order ters. of the scattered particle 
wave function, with an arbitrary scattering pote: .ai. 


Introductica 


This paper deals with two methods of :reating scattering problems in 
quantum mechanics which differ widely in approach. One, the well known 
Born approximation, is a perturbation theory solution of the Schrédinger 
equation. The other, described in recent papers by R. P. Feynman (1, 2), 
deals directly with the solutions of the wave equation, rather than with the 
equation itself. 


In view of the fact that Feynman's approach leads to a simplified formu- 
lation of quantum electrodynamical problems, it seems desirable to establish 
a link between this procedure and the familiar Born approximation. The 
present paper accordingly carries out explicitly Feynman’s treatment of a 
simple scattering problem and demonstrates the equivalence of the two 
methods. 

The Born Approximation 


Suppose that a particle of mass m and initial momentum hk* is scattered by 
a potential field V(7). To apply the Born approximation, we set v(7,%) = 
u(7)e““, with w = hk?/2m, and proceed to solve the time-independent 
Schrédinger equation 


{v2 + k? — (2m/h*) V (7)} u (7) = 0. 


The method is well known, and we include here only those equations which 
will be compared directly with the Feynman solution. The wave function for 


the incident particle is u = e**-?. Green’s solution of the first order equation is 


(72) = (2m/h*) f G(Fs, 71) V (Fi)uo(71)dr1, (1) 
where the Green’s function G(72, 71) = — (4% | 72 —71 |) etlrs— Tl satisfies the 
eccnicies {vst + } G(Fs, 71) = 8(F2 — 71). (2) 


1 Manuscript received August 8, 1950. ‘ 
Contribution from the Department of Physics, University of British Columbia, Vancouver, 
B.C. This paper is to be submitted in partial fulfillment of the requirements for an M.A. degree. 
The research was carried out with the aid of a Scholarship from the National Research Council 
of Canada. 
*The symbol h represents Planck's constant divided by 2m (h = 1.054 X 10-*" erg. sec.). 
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We write d? = dxdydz; unless otherwise specified, the integration is to be 
taken over all space. The subscript 2 on V? means that it operates on the 
variables 72, and 5(72 — 7:) = 6(x2 — x1) (v2 — yu) 5(z2 — 21). 
For solutions at large rz we make the approximations 
| 72 — Fi | = 42 — Fi- Fo/r2 = re (3) 
in the exponent and denominator respectively of the Green's function. The 
first and second order solutions are then: 











ikr we 2 ‘a 
u(%2) = — a "| vr nett-to Far, (4) 
2rh? T2 
2 pikr > +. ght ke ee Ri rn 
u2(72) = ( “ ) a | veri? oo | vend Far, dr; (5) 
2rh? 12 3 


where ky = k¥2/r2, kg = RP s/rs. 


Feynman’s Formulation 


Feynman’s method is based on the probability amplitude for a particle to 
move from one space-time point to another. This function, called the kernel, 
is a sum of contributions e*5/* from each possible path, where S is the classical 
action along the path. Denoting the point 7;, fi: by I and 72, te by 2, K(2, 1) 
is the amplitude for a particle known to be at 71, at the time ¢; to arrive at 72 
at the time fz. The wave function ¥(2), representing the amplitude for the 
particle to be found at 72, tz, is then given by the expression 


¥(2) = f K(2,1) ¥(1) dri. (6) 
Feynman shows that the kernel is that solution of 
{ihd/dt, — H2} K(2,1) = ih8(2, 1) (7) 


which vanishes for fz < t;, where H is the Hamiltonian operator for the particle 
(the subscript 2 indicating that it acts on the variables 2), and 6(2,1) = 
5(7%2 — 71) 5(t2 — t:). The kernel for a free particle is shown to be (for t2 > 41) 


= m 3 im \F2 = Li 3 8 
0 lac * 5 exp| (h — h) @) 


For a perturbation theory treatment of a particle in a potential V(7), 
Feynman shows that the kernel may be expanded in increasing powers of V: 


K(2, 1) = Ko(2, 1) + Ki(2, 1) + K2(2,1) +..., (9) 
where K,(2,1) = — (é/h) f Ko(2,3) V(3) Ko(3, 1) drs (10a) 
K;(2, 1) = (—1/h)? ff Ko(2, 3) V(3) Ko(3, 4) V(4) Ko(4, 1) drs drs (100) 


etc., using dr =drdt. The time integration is limited to the interval 
te > ts > tr in (10a), te > ts > t4 > th in (105), by the condition K(2,1) = 0 
for tle < ky. 
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The decisive advantage of Feynman’s theory is the interpretation of the 
above formulas in terms of successive scattering. In (10a) for example we 
view the particle as moving freely from 1 to 3 (amplitude Ko(3, 1) ), being 
scattered by the potential at 3 (amplitude —(#/h) V(3) dr3), and moving as a 
free particle again to 2 (Ko(2,3) ). This is summed over all points 3, thus 
giving the term K,(2, 1). The successive terms of (9) are regarded as ampli- 
tudes for the particle to be scattered 0, 1, 2, etc., times by the potential V. 


Corresponding to the expansion (9) of the kernel, the perturbed wave 
function may be written 


¥ (2) = Yo(2) + vil2) + ¥2(2) +..., (11) 

where yo(1) is the wave function of the incident free particle, and 
¥o(2) = f Ko(2, 1) yo(1) dri, (12a) 
vi(2) = —(i/h) J Ko(2,1) V1) Wo(1) dri, (126) 


etc. 


Comparison of the Methods 


From the foregoing discussion it is possible to show a close correspondence 
between the two methods. The similarity of Equations (12b) and (1) suggests 
identifying ' 

G(72, File? = —(ih/2m) J Ko(2, Ie" dh; . (13) 
Operating on (13) with {ihd/dt2 + (h?/2m) V2*} gives on the left, by (2), 

(h?/2m) {V2 + k?} G(F2, File *** = (h?/2m) 6(72 — Tie", 
and on the right, using (7) with H = — (h?/2m)V,’, 
—(ih/2m) f ih 5(2, 1)e*" dty = (h?/2m) 6(72 — 7i)e**". 

The equivalence of these two expressions does not prove Equation (13), since 
an arbitrary solution of the Schrédinger equation can still be added to one side. 
It does show however that there is an exact parallel between the two methods 
in spite of the wide difference in approach. Without giving a full proof of 


Equation (13), we will study the correspondence of the two methods from 
another point of view. 


The n** order perturbed solution is obtained by using for the kernel the first 
n + 1 terms of the series (9). We now seek to modify the Hamiltonian of 
the particle by introducing a potential function V,(7) such that Ko(2, 1) +... 
+ K,(2, 1) is the exact kernel for a particle moving in the potential V,(7). 
Applying Equation (7) gives 
{ihd/dt, + (h?/2m) V22 — Vn(2)} {Ko(2,1) +... + Kn(2,1)} = ah6(2, 1), (14) 
which defines the function V,(7). 

A general expression for the n™ order kernel is 


K,(2,1) = —(t/h) f K,(2, 3) V(3) Kn-a(3, 1) drs. 
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Using the free particle form of Equation (7) gives therefore 
{ihd/dt, + (h?/2m)V.2} Kp(2, 1) = V(2) Kn-1(2, 1). 
Thus we see that I’,(2) is given by the equation 
Vn(2) {Ko(2,1) +... + Kn(2,1)} = V(2) {Ko(2,1) +... + Kn-a(2, I}. 
Multiplying this equation by o(1) and integrating over d71, we obtain 
Vn(2) {Yo(2) +... + vn(2)} = V2) {Yo(2) +... + vn1(2)}. (15) 
From the definition of V,, the wave function yo +... + Yn must satisfy 
the modified Schrédinger equation containing this potential: 
{ihd/at + (h?/2m) V2 — Vi} {yo t+... +n} = 0 
Using the expression (15) for V,, this becomes 
{iha/at + (h2/2m) V2} {yo t... +n} = Vidot... + ona}, 
and subtracting the corresponding equation with m lowered by one gives 
{ihd/dt + (h?/2m) V?} wn = V War. (16) 
This is exactly the n order equation of the Born approximation. We have 


thus proved that the Born and Feynman solutions are identical, provided that 
the same incident particle wave function pp» is chosen. 


The condition for validity of either approximation is that the terms of the 
series expansion of the wave function (11) should decrease rapidly, i.e., that 
|¥1|<«]|yo| = 1. Since the two methods give the same result for yj, it is 
clear that they are valid under the same conditions. 


Explicit Solution by Feynman’s Method 
To illustrate the equivalence of the two methods, we will use Feynman’s 
formulas to calculate explicity the first and second order terms of the wave 


function. We set Yo (1) = e#-7:— 4) for the wave function of the incident 
particle. The zero-order wave function (12a) should of course be the same. 
This is easily verified by direct integration using the expression (8) for the 
kernel. 


The first order term, from (120), is: 
‘Tr ! m ; |= \7s ee 4 see 
2 eat ead 2X See eee ee et(k-ri—wh)g dt 
es fe Hie l ka e ss 


= —(i/h)eilsf T(|%_ — 71|) VF de® Fd7 





where T(r) -| (Qriht/m)~2eimr?/2htgitdy 
0 


(usingr = | 72 ” Ti | ,¢ = tg — t;) is the integral on the right of Equation (13), 
multiplied by e'!:, 
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Now according to Feynman’s interpretation of the perturbation formulas, 
¥i(2) is due to particles scattered once at 1 and then proceeding as free particles 
to 2. Since a free particle has a definite velocity hk/m, we should expect a 
contribution to T(r) only from the neighborhood of the point t = (m/hk)r 
We rearrange the exponent, and make the substitution ¢ = (m/hk) r (1 + £): 


T(r) = (=.)'[r texp Ea — +} br = : ‘la 
2xth 


= ia) GY [1p + oto (+3 +8) Ia 


a. = (zye |" (1+ gy-lebibe/A1+ Bar. 

2Qrih\2xi/ +r 

The coefficient of zt in the exponent is $kr¢/(1 + &). For fairly energetic 
particles at large distances re, kr > 1. Hence except for the region § < 1 the 
exponential is a very rapidly oscillating function of — The function (1 + &)~? 
varies slowly (except at — — 1, and here the exponential oscillates with 
infinite frequency), and hence the integral vanishes except for very small 
values of ¢. Thus the physical reasoning of the previous paragraph is con- 
firmed. For small , we take 1 + ~ = 1, and then, since kr > 1, extend the 
lower limit to — ©. The integral becomes* 


T(r) - = (£) <I —(kr /2i) Fe 
2xth \2xi/ or J-o 


= (m/2nih)e**/r. 
Comparing this with the Green’s function, we see that Equation (13) is verified. 


Applying the asymptotic approximations (3), the first order wave function 


becomes og re 


wt ) z an _ > 
¥i(2) = — — tet [Viineid- in 77, (17) 
2rh? TP? 
which agrees exactly with the Born approximation result (4). 
From (10d), the second order wave function is: 


¥2(2) = (—i/h)* ff Ko(2,3) V(8) Ko(3, 1) V(1) Yo(1)drs dri 
= (—i/h)? ef f T(|7.—75]) ao T (|%:—73|) VF) ef -FidF a7, 


i(kr2—wts) 
“a (# ay eters % = [ve sen tee: 73£ a lve ) ilk hs): id? dr; (18) 
t To 3 








using again the approximations (3). This is identical with the Born result (5). 


* In spite of the approximations used, this result is exact. Withx = 1+, p = —ia = 
— hikr, the integral is ta 3/2e—px+ia/xdx, This can be evaluated by using a table of La- 
place Transforms, such as (8). The method used has the advantage of affording a physical inter pre- 
tation. 
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It will be noted that the successive integrations involved in the higher 
order approximation do not lead to a more complicated time integral (as might 
be expected), but rather to simple products of the same function T(r). This 
is due to the separation of the time dependence of Equation (13) in the factor 
ez, It is clear, in fact, that this time integral of the kernel will always give 
the same function provided only that the wave function can be separated in 
this way. Thus this integration, the most difficult step in the Feynman 
calculation, need not be repeated in every problem. 


Acknowledgments 


The author wishes to thank Dr. F. A. Kaempffer for suggesting the topic 
of this paper and for helpful discussion, and also the National Research Council 
of Canada for the award of a Bursary under which the research was carried out. 


References 


1. FEYNMAN, R. P. Rev. Modern Phys. 20: 367. 1948. 

2. FEYNMAN, R. P. Phys. Rev. 76: 749. 1949. 

3. MacGnus, W. and OBERHETTINGER, F. Special functions of mathematical physics, Chelsea 
Publishing Company, New York. 1949. p. 127. 








72 


UPPER AND LOWER BOUNDS FOR THE ASYMPTOTIC NEUTRON 
DENSITY IN MILNE’S PROBLEM FOR THE SPHERE! 


By W. SEIDEL? AND R. E. MARSHAK? 


Abstract 


Rigorous upper and lower bounds for the asymptotic neutron density in 
Milne’s problem for the sphere are calculated by means of an extension of Hopf’s 


method for the plane. 


We consider the following problem: a black sphere of radius a, i.e., a sphere 
which absorbs completely all neutrons which fall upon it, is surrounded by an 
infinite noncapturing medium which scatters neutrons isotropically without 
changing their velocity. No sources are present in the outer medium except 


; F . 7 io 
that a current density of ae (F is a constant) is assumed to exist in 
™r 


the direction —r. We wish to determine upper and lower bounds for the asymp- 
totic neutron distribution in the (outer) medium. 


An exact determination of the neutron density for the present problem—by 
transform methods or otherwise—seems out of the question (Berlin (1) and 
Hansson (2)). Fortunately, for some problems of practical interest it is 
sufficient to know the asymptotic behavior of the neutron density. In a pre- 
vious paper (Marshak (4)), * one of the authors solved the transport equation 
approximately by expanding the angular distribution in a finite number of 
spherical harmonics; this method has the advantage that it yields an approxi- 
mate expression for the neutron density over the entire range. The disadvantage 
of the spherical harmonic method lies in the difficulty of estimating the error 
of any particular approximation. To check this point for the sphere, we derived 
the equivalent integral equation and from it — using methods first developed 
by Hopf (3)—we now derive rigorous upper and lower bounds for the asymp- 
totic neutron density. A comparison of the results with the results of the 
spherical harmonic method enables one to judge the accuracy of the successive 
approximations in the harmonic method. A further value in reporting the 
present calculations lies in the possible extension of this method to other 


problems. 


As shown in I, the asymptotic solution for [rpo(r)] is of the form B(r — ro), 
where ro is a function of a and for large a should go over into (a —0.7104). We 
start with the integral equationf for [ryo(r)], which is equivalent to the tran- 


1 Manuscript received November 7, 1949. 
The contents of this paper were contained in a Montreal Report (MT-49) dated April 15, 


1944. 
2 University of Rochester, N.Y., U.S.A. 
* This paper will be referred to as I and the same notation will be used. 
t The derivation of this integral equation is straightforward and is given in our original report 


MT-49; cf. also references in (1). 
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sport equation plus boundary condition y(a, u)=0 for u >0, and from this 
integral equation we obtain rigorous upper and lower bounds for ro. 


The integral equation in question* is: 
co siti a 
ryo(r) = 5 | r'po(r’) {E(| r—r|)-E (Vr —-@+Vr? — a?) } dr’. (1) 


To derive upper and lower bounds for the asymptotic neutron density, we 
write: rfo(r) = r + q(r), where g(r) satisfies the integral equation: 


ore z|. g(t){E(|r — tl) — E(Wr — a? + VP — a4) at ee , (2) 
with E3(r) = — {aE.(r — a)— Ex(r — a) + E; (Vr? — a’) }. (2a) 


Making use of the properties of the E,-functions, it can be shown that E;(r) 
is negative for all r > a (cf. Appendix). Further we remark that g(r) satisfies 
the inequality —a < g(r)< O(cf. Appendix) and must be a monotonicallyt 
increasing function of r. We now integrate Equation (2) with respect to r from 
r2atoM>r. We get: 


M 1 M ‘co 1 M co 
| q(r)dr = il ar| q(t)E(\r — t|) dt — 2/ ar| q(t)E(V ra? 


conasilatee M 
+V 2 — a!) at + + | Ex(r)ar. 
In the second and third integrals the order of the integration may be reversed 


since the integrals with respect to ¢ converge uniformly for r in r& r< M. 
Thus we find: 


M 17 M 1 [¥ ra 
| q(r)dr = 2| athat | E(r —t)dr + 3] q(e)dt | E(r — t)dr 
1 M t 1 oO M 
+—| q(t)dt| E(t—r)dr+ i/ a(t | E(t — r)dr 
2 rT r 3 M r 
1 co M etisalat 1 M_ 
= il q(t)dt f E (Vr —@@tvPf=— a’) dr + | E;(r)dr. 


With the help of relations between £,-functions, we obtain: 


* The function En(x) ={> é **u—”" du for x > 0; the E,-function is usually written without a 
subscript. For the elementary properties of the functions En(x) the reader is referred to E. Hopf, 
loc. cit. pp. 26-27. : ’ 

t This result follows from physical considerations and can be seen very roughly as follows: if 
the black core had a smaller radius, the asymptotic solution for Wo(r), namely( 1 —"°\, would be 

r 


valid further inward than tt is actually valid. Since this is not the case, the density must drop more 
rapidly than om by ro because neutrons which further outward are scattered back are now—in 
the neighborhood of the core—lost owing to the absorption by the black core; this effect becomes more 


and more important (monotonicity!) as the core is approached. We have not f er, 
an analytic proof of the monotonicity. - , succeeded in obtaining 
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M Att 1 M 
| q(r)dr = 4 [ow { Eo(r — t) — Ex(M —1t)}dt+ > | g(t) {1 — Ex(M —#) }dt 


M ‘<0 
+ Lf q(t) {1 — Ex(t -aja+t/ g(t){ Ex(t — M) — E2(t — r)} de 
2 Jr 2JmM 


© M M 
— il a(t [ EVr-a&+Ve- a®) dr +4] E;(r)dr. 


Hence, os 
0= [aw { E2(r —t)— E,(M — t)} dt -| q(t)E2(M — t)dt 


M oo 
-| g(t)E2(t — r)dt +] a {E.(t — M) — Ex(t — 1)} dt 


- [awae [ee r? —@+Vf—a’)dr+ | Baar (3) 


There are eight integrals which will be denoted by J;, I2,..., Is respectively. 


We now examine each integral as M > o.* It is clear that J; remains un- 


co ‘co co 
changed, J2 0, 14-3 — | a(t) Ex(t —1)dt,I,—0, I; -| attat| E (Vr?—a? 


+Vp— a’) dr, h| E;(r)dr. It remains to study J; and Js. 


v 


We shall now use the fact that |q(t)| is a monotonically decreasing function. 
We then have, since \a(t)| is bounded, that 


lim | g(t)| = ro 
t—-@ 


exists. We consider: 


© M 
lim | | q(t)Ex(t — M)adt -| g(t)Ex(M — pat | 
M T 


M-—© 


“= ! | “or nee = [ace , »)Ex(x) | 
M-+@ 0 0 


me 


© M 
=— rol E2(x)dx — lim | q(M — x)E2(x)dx. 


We wish to show that: 

M-r 
lim | Q(M — x)E,(x)dx = — (4) 
0 


M-© 


To this end, we write: 





M-r 


|q(M — x)|Ex(x)dx > Fo | Bale) - nf ~ Eat = at. 
0 





* The passage to the limit under the integral sign may be readily justified in each case. 








I 
o 
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On the other hand,* 
M-r 
| q(M — x)|E2(x)dx 


M M/2 M 
= [ | g(t)| Ex(M — t)dt = I. | q(t) | E2(M — t) dt + {ha | Ex(M — t) dt 


<a(Dleo|(4-»)* |e) fa 
<o)o(4—»)+ [DIE 2 


Thus, we have 


na — E,(M — »} < [lacar — x)| Ex(x)dx < Ey (*)0(¥- r) 


+leSite- =G)} 


Allowing M — in the above inequality, we obtain 














M-r 
< |im \, \q(M — x)|E2(x)dx < —. 
2 M—+© . 
This proves (4) so that finally 


lim (J3 + Js) = 
M-@ 


Letting M — = in (3), we obtain: 


én | a(OEx(s hile [emese oe [ea [2 (VA — a? 
+V# — a’) dr+ | “Exar (5) 

Let + = a in (5); then: 

0 =| ones — a)dt +] acnae[ BO/ rr —at+VP—a)dr - | Burra (6) 

If we multiply (6) through by —1, we find 

| “la as{ “EVA a@+Ve—at)dr + | “lao E,(t — a)dt = [zslar, 

Since a monotonically decreasing, we may write: 


ro} | ” Balt —a)dt + | “dl EW r— a? +V# — a’) ar\ < [een dr, 


* We implicitly assume T < — this is justified since we allow M — © in the end. 
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or finally: 
ad — 
[| Buel a 
Bl, 
1 — 
i+] at| E(V ? — a + V8 — a?) dr 
@ 


a 


(7) 


The inequality (7) gives a bound from above for the quantity ro. To obtain a 
bound from below, we return to Equation (5), replace there + by u, and inte- 
grate with respect to u from 7 to M: 


M 1“ M © 
0 -| au| q(t)E2(u — t)dt -| au | q(r)E.(t — u) dt 
- © oo iit eaten, M oO 
-| au | ata | E(Vr — a + V#— a’) dr +] in| E;(r) dr. 
Reversing the order of integration, as before, 
r M M M M ¢ 
0 -| athat| E,(u — t)du +] attat| E.(u — t)du -| athat| E.(t—u)du 
‘co M M co oc 
-| a(bat| E2(t—u)du -| au ata] Er —a? +V2 — a’) dr 
M r T a u 


M co 
a | au | E;(r)dr. 


Again using identities on E,-functions: 


r M ( 1 
0 -| g(t) {Eals =f) bia ok dt +| q(t) a ~E(M— ok dt 
-| a {3 te ap dt -| {Est «MR aha (8) 
-| auf acnatl ev r—ar+ Ve— a®) dr + [ au] Barna, 


Allowing M — as before, the only integrals that require special study are: 


lim | - | “q(t)Ea(M — dt - | aneat — M) at] 
M0 tr M 


= lim Bra — u)E;(u)du ~ | acar + wExuau |. 


M-—o 0 


By an argument similar to that used in deriving (4) we obtain the limit + > To. 


Hence, on passing to the limit M — o in (8), we get: 


= | ‘q(t)Es(r — t)dt + [ ames — r)dt + = 


-| au [a (oat| "E(V e+ VF — a*)dr + i | Baar, (9) 
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or setting r = a, 


0 =| aes —a)dti+ =r -| a [ awa [Be rP—g@4+Vef— a*)dt 


“fb | au| E;(r)dr. 
Multiplying by —1, we obtain: 
=n = | \q(t)|Es(t — a)dt -| lawl [ au| E(Vr — a +V#— a?)dr 


F [a [Tze (10) 


Now, 
| E(Vr —a’ +VP —a?)dr <| E(r—a+V2—a?®)dr = E,(u—a+Vf—a?) 


and 


[au E(Vr -v@+VP- a)dr < E,Ve —a’) < E;(t — a). 
Hence, we find from (10): 
+ ro 270 \; _ | in| du| E(Vr -a+ Ve a?) ar + dul |E3(r)|dr, 


wt 


and finally: 
| au | |E,(r)|dr 


Pg ea anh aetna 2 (11) 
a+ al au] EW r —-a+V 2 — a) dr 


From (6) and (10) one may obtain upper bounds for 79 which are better than 
(7). We have, indeed: 


Faw E.(t — a)dt = | [ea dr — J iacola EG —a?+V# —a®)dr, 


and 
| \g(t)| E(t — a)dt = = +] iacolar| aul EVP —-a@+V2—a)dr 
_ | au| |E3(r)|dr. 





This, in turn, yields: 


, eee , Oe  caeetee acacasiel 
| \q(t)| Eo(t — a)dt < | \E3(r)|\dr — ra aul EV? — a +VP —a’*)dr, (12) 
a a a a 


and 
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[locolesc —a)dt > {2 + [ae] dul evra et + Va—a)arh re 
-[ au [ [Bacar (13) 


Since E.(t — a) > E;(t — a), on subtracting the first equation above from the 
second, we find: 


io {Ex(t —a) — Ex(t — o)}dt>ro} > 4 [a | eeame are. 
+ [a au ee/e=o + Vea ar} 


- [ lzsnar [au | TBscoler, 
and finally: 
| [esa + dul |E3(r)|dr 
SSeS Oe’ 
2. + ar] E(V r—a?+ V?—a*)dr +| a] au E(V ra? + V?—a*)dr 


(14) 
The second method is based on the fact that: 


E,(t —a) -— = Ext —a)>0, 


and that | q(t)| is‘monotonically decreasing. We find that: 


[, {2 —a)- = Balt a) aco) <0. 


a 


Integrating by parts and observing that the integrated part vanishes at both 
limits, we obtain: 


| [lac 2 E,(t — a) — Ex(t — ab ae > 0. 
Using (12) and (13) again, we obtain: 
0< = [Leslee + 3 ofa [Berm + Vf — a?) dr 
“f+ [fe =e ov} [el 


and finally: 
foe} co ce 
; | |E3(r)|dr +{ du |E3(r)!dr 


6 Sg a anne 
2! as| ECVP a4 VE—aildr + 2 + [tel dul BO Aa VA —alar 
(15) 
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It is to be expected—and has also been checked—in analogy with the plane 
case that (15) is a more refined upper limit than (14), which itself is a more 
refined limit than (7). We have tried to obtain a more refined lower limit for ro 
than is given by (11) but without success. It seems impossible to do so in 
principle since the part of the integral kernel which depends on the square root 
always makes the inequalities go in the wrong direction. Since more delicate 
properties of the difference of two £,-functions are employed in the derivation 
of the refined upper limit (15) than in the derivation of the crude lower limit 
given by (11), it is clear that the value of ro given by (15) must be closer to the 
actual value than the one given by (11). This is made very plausible by con- 
sidering the limit a — ; the integrals involving the E,-functions of square 
roots then ey zero and the bounds for (ro — a) — it is more revealing 
to talk about (ro—a) than r>— are given simply by: 

3 17 

= < (ro — a) < - (16) 
The correct value of (ro—a) is —0.7104; it is seen that the upper limit is out 
only by 0.3%, while the lower limit is in error by 5.5%. Since all the integrals 
involved are monotonic functions of a, it is reasonable to suppose that a similar 
situation prevails over the entire range, except for sufficiently small* a 


The integrations in (15) involving the square roots can only be performed 
by numerical means. It is convenient to rewrite (15) in the form: 


17 2 2A 
“nota c)e(tt0) 
a a(2 + + - 


(ro — a)y 2 ’ (17) 
— +(2 B+ c) 
3 
where 
oc 
A -|. E;(V# — a)adt, 
B =| aul EWr — a + V# —a®)dr, 
¢ -[ in| au 2 EVe-a@ —-a+ Vr —a? — a*)dr, 
D -| aul 2 E;(V ee a*)dt. 
In terms of the above notation the lower bound for (79— a) is: 
_ eS —aC+D 
(ro — a), > —— pomeneere (18) 
—+C 
* For very small a, Equation (15) predicts an a? log a behavior for roas a — O, whereas the correct 


behavior is a®?. C “uriously enough, Equation (11) gives the correct behavior of ro as a + O. 
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The integrals A, B, C, and D have been computed for a = 0.5, 1, 2, 5 and are 
given in columns 2-5 of Table I. Columns 6 and 7 contain (r9>— a), and (ro— a); 
respectively. 


TABLE I 
a A B Cc D (ro— @)u (ro— a)i 
0.5 0. 2071 0. 1764 0. 1180 0. 1452 —0.340 —0. 363 
1.0 0. 1516 0. 0932 0. 0560 0.0941 —0. 503, —0. 544, 
2.0 0. 0968 0. 0377 0.0174 0. 0457 —0.630 —0. 681 
5.0 0.0442 0.00782 0.00187 0.0120 —0. 692 —0.740 
© 0 0 0 0 —0. 708 —0. 750 
Appendix - 


We first prove a lemma: 

Lemma. For all r > a: 
Ex(r) = — [aE.(r — a) — Ex(r — a) + Es (V7? — @)]< 0. 
By the mean value theorem: 
E;(r —a) — E; (WV? —a?) = E',(¢) [r —a —- V7? —'] = E(t) [V7?— a? — (7-2), 
where r-a<t< V7 — a. 
Now, 
aE2(r —a) — Ex(r —a) + Es (V7? — a?) =a E2(r — a) — Ea(t) [W7? — a? — (7 —2)] 
> aE2(t) — E(t) [(W7? — a? — (r — a)] = Ex(8) [r-W7? — a] >0. 

This proves the lemma. 


We next prove that if f(r) is the solution of the integral equation: 


f(r) = 5 [70 {E(|r —¢]) —-EWr—a + V2 —a)} de, 


so normalized that f(r) = 7 — a + g(r), then 0 < 9(r) <a for allt > a > O. 


We first observe that since: 


|r7—t|< V2—a?+V#— a? and E(t) is a monotonically decreasing function, 
E(\r — |) - E(V?— a2 +V# — a*)>0. Next, evaluating the integral: 





4 ede -— 1) -EVP-@4+Ve-a')}dt =1 — 5 Bale — a) 
2 - | EW A= a+ Via") dt 


we obtain: 
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oo 
1 -i/ {E(|r ~ t|) = Ew r -—a+ V P =a!) | dt + = Ealr — 0) 
oo ee es 
+ i/ E(V? —a+V# — a’*)di. 


Secondly, evaluating the integral: 
| (t—a){B(|r — |) — EWP oat + VP mo) at, 
we obtain: 


pean 1 lu afl ~ |) —-EWenat +VP— a8) Jat 


co Siisinaieieeceliaeens 
+ 1 (t -—a)E(V 7 -a@ +VF a*)dt — 5 Ear — a). 
Furthermore, 
@ _— a 
1} tE (V7 —a’ +vVP — a*) dj = a (V7 By a’), 
so that the above equation becomes: 


(2) r-a= 3 { —a){E(|r —t|) —-EVP2-a&+V# —a*)} dt 


~ = Eals —a)+ ~ Ea —a’) — 2 | EW? =~ gv pa a*)dt. 


Multiplying Equation (1) by a and adding to (2), we get: 
(3) r= + [i tE(e -—t|) -EW?P-@+Ve—a*)} de 
+> (aEs(r — a) — Ex(r — a) + Ex(V7? —a)}. 

Since in (2): 

E;(r — a) — E; (V7? — a?) + a E(V?—a?+V# — a*)dt>0, 
we obtain: 
(4) +r-aK< Lf (t—a){E(|r -—¢]) -EWP—-a@&4+V2—a*)} de 
and, in view of the lemma, from (3) we find 


(5) r> 4 [ite(r -A) —~E(V? — a2 +V 2 —a?)} dt. 


We now introduce the sequence of functions: 


(6) filr) = 7 —a,fn(r) = 4 hau {E (|r — tl) 
~E(VP-a+V# —a’)} dt. 








82 CANADIAN JOURNAL OF PHYSICS. VOL. 29. 
Inequality (4) shows that: 

fi(t) < fa(r). 
Multiplying both sides by = {E(|r—¢) -EWP-a@ + Ve —a)} and 


integrating, since the multiplying factor is positive, we find: 
f2(r) < fa(r) 
and generally: 
(7) fn-a'7) < fa(r)- 
Returning now to inequality (5), from +r — a < 7 and the positivity of the 
kernel, we find: 


f2(r) <1 [ifets _ t|) - E(v ?—a?+V2—a?*)} dt <7. 


Repeating the argument, 


fa(r) -1 |i {E(|r —t)- EW? -—-?+v? — a’)} dt 
<1 [ited — 1) -EWer-a@+Ve—a)} dt <r. 


Thus, 
(8) fr(t) <r. 
In any finite interval, therefore, the sequence { f,(r)} is bounded and ac- 
cording to (7), the sequence is monotonically increasing. Consequently, the 
limit: 

lim f(r) = f(r) 

"m—» 0 
exists. By well known theorems we are justified in taking the limit n > « 
under the integral sign in (6), which shows that f(r) is a solution of the integral 
equation: 


f(r) = 1 lho {E(|r —¢) — E(V? —-a+vVePr — a*)} dt. 


It follows, moreover, from (7) and (8) that the solution f(r) satisfies the 


inequalities: 
r—a<f(r) €r. 


Applied to g(r), this means that 0 < g(r) < a, as was to be proved. 
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A METHOD OF MEASURING GAMMA-RAY ABSORPTION 
COEFFICIENTS AT 0.51 Mev! 


By P. E. ARGYLE, G. M. GRIFFITHS, AND J. B. WARREN 


Abstract 


The coincidence method described makes use of the angular correlation 
between pairs of annihilation quanta to eliminate the effect of scattered radiation 
usually present in absorption measurements. Mass absorption coefficients in 
(gm./cm.?)—! found by this method are as follows: 

Lead 0.1519 + 0.0008 
Copper 0.0819 + 0.0005 
Aluminum 0. 0831 + 0.0013. 


Introduction 





In the usual methods (1) employed for accurate measurements of gamma-ray 
absorption coefficients, the main difficulty is to prevent scattered radiation 
from reaching the detector. We have circumvented this difficulty by using 
the precise 180° angular correlation (2) between two annihilation quanta to 


measure the absorption of 0.51 Mev. radiation. 


Two gamma-ray scintillation counters were placed in line several meters 
apart, one on each side of a thick source of positrons. The coincidence rate 
between the counters was measured as a function of the thickness of the 
absorber which was placed between the source and one counter, as close to the 


source as possible. 


Anthracene Anthracene 
Crystal Crystal 
ge sae 





Coincidence 
Output 
Fic. 1. Apparatus for measuring absorption coefficients at 0.51 Mev. 


Counter A (Fig. 1) acts as a collimator for the gamma-ray ‘beam’ from 
the source to counter B, but does not itself scatter a significant amount of 
radiation back into B. Radiation scattered by the absorber and by the 
laboratory walls, while it may enter one of the detectors, cannot produce a 
coincidence count in both detectors owing to the exact antiparallelism of the 
pair of annihilation quanta. This is a simple method of producing a beam of 
monochromatic gamma-rays collimated to any desired degree of accuracy— 


1 Manuscript received August 28, 1950. 
Contribution from the Department of Physics, University of British Columbia, Van- 


couver, B.C. 
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which could be used, for example, to measure the variation in efficiency of a 
Geiger counter over its volume. 


Experimental 


The scintillation counters consisted of blocks of anthracene placed on 
RCA 5819 photomultipliers, and the source of annihilation radiation was Cu 
produced in the Chalk River pile. The distance between each detector and 
the source was 220 cm. and the plane angle subtended by the phosphors at the 
source was just 3/4°. The pulses, after amplification, were standardized by 
fast multivibrators with dead times of 17 usec. and then fed into a Rossi circuit 
of resolving time 0.15 usec. Considerable care was taken to keep the over-all 
gain of the system, and hence the efficiency of the counters, constant during the 
runs. The stability was checked by taking readings on the coincidence rate 
with no absorber and comparing the half life deduced from these readings 
with the accepted value. With the multiplier operating conditions chosen 
the noise rates were about 10 per second, which gave an insignificant number 
of coincidence background counts. 


Results 


The experimental readings were corrected for accidentals, dead time losses, 
and source decay. The value of half life used for the decay correction was 
12.88 + 0.03 hr. (5). The table below shows the experimental absorption 
coefficients obtained from our results by the method of least squares, the 
values predicted by theory, and values obtained by other experimenters. The 
photoelectric absorption coefficients were obtained from the curves of Hulme 
and others (4), and the Compton absorption was calculated on the basis of the 
Klein-Nishina formula. The agreement with theory is within the errors 
quoted for the theory and the statistical errors of the experiment. All results 
are given in units of (gm./cm.?)7!. 


: p | Total a 7 
Absorber | Photoelectric] Compton | theoretical | Present authors | Glendenin} Cowan 





























absorption | absorption) absorption Ref. 3) | (Ref. 1) 
Pb 0.0802 | 0.0675 0.1477 | 0.1519 + 0.0008] 0.154 | — 
Cu | 0.0028 | 0.0778 | 0.0806 | 0.0819 + 0.0005; -— | 0.079 
AL | 0.0002 | 0.0823 } 0.0825 | 0.0831 + 0.0013 | 0.082 | 0.083 
' | | i 





The maximum thicknesses of absorber used were 4 cm. for lead, 4.7 cm. for 
copper, and 5 cm. for aluminum. Fig. 2 shows the counting rates in the two 
counters and in the coincidence channel, after corrections, for the case when 
lead absorbers were used. 

These curves indicate clearly how the effect of scattering from walls, tables, 


and supports distorts the counting rate curve for counter B when the absorber 
is put between that counter and the source, while the coincidence curve 
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remains straight. The stability of the apparatus is indicated by the fact that 
the counting rate in counter A remained constant to within 3%. The actual 
stability is probably better than this, since a considerable fraction of the 3% 
rise in counting rate of ‘A’ as the absorber thickness is increased can be attrib- 
uted to backscattering from the absorber. 


yy canal 
fe 





counter B 


Log Counting Rate 


Absorber Thickness 


om,/cm* 


Fic. 2. Channel counting rates as a function of absorber thickness. 


The aluminum used was Alcan 2S, which the manufacturers state is about 
99.5% pure. The copper was made electrolytically. The purity of the lead 
was better than 99.9% 
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